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ABSTRACT 

We have searched for in-faUing steUar streams on to the local Milky Way disc in the 
CORAVEL and RAVE surveys. The CORAVEL survey consists of local dwarf stars 
(Nordstrom et al. Geneva-Copenhagen survey) and local Famaey et al. giant stars. We 
select RAVE stars with radial velocities that are sensitive to the Galactic vertical space 
velocity (Galactic latitude b < —45°). Kuiper statistics have been employed to test 
the symmetry of the Galactic vertical velocity distribution functions in these samples 
for evidence of a net vertical flow that could be associated with a (tidal?) stream of 
stars with vertically coherent kinematics. In contrast to the 'Field of Streams' found in 
the outer halo, we find that the local volumes of the solar neighbourhood sampled by 
the CORAVEL dwarfs (complete within ~3 xlO"" kpc^), CORAVEL giants (complete 
within ~5 xlO~^ kpc'^) and RAVE (5-15 per cent complete within ~8 kpc'^) are devoid 
of any vertically coherent streams containing hundreds of stars. This is sufficiently 
sensitive to allow our RAVE sample to rule out the passing of the tidal stream of 
the disrupting Sagittarius (Sgr) dwarf galaxy through the solar neighbourhood. This 
agrees with the most recent determination of its orbit and dissociates it from the 
Helmi et al. halo stream. Our constraints on the absence of the Sgr stream near the 
Sun could prove a useful tool for discriminating between Galactic potential models. 
The lack of a net vertical flow through the solar neighbourhood in the CORAVEL 
giants and RAVE samples argues against the Virgo overdensity crossing the disc near 
the Sun. There are no vertical streams in the CORAVEL giants and RAVE samples 
with stellar densities >1.6 xlO"' and 1.5 xlO'^ stars kpc^'^ respectively and therefore 
no evidence for locally enhanced dark matter. 

Key words: solar neighbourhood; Galaxy: disc; Galaxy: formation; Galaxy: kine- 
matics and dynamics; Galaxy: structure; dark matter 
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I INTRODUCTION 

Tidal streams are filaments of debris, containing stars, gas 
and possibly dark matter (DM), that have become stripped 
from a disrupting satellite dwarf galaxy (or star cluster) 
along its orbit with respect to the centre-of-mass of a more 
massive galaxy. The instantaneous position and velocity of 
a star and the Galactic gravitational potential determines 
its past orbit. For stars in tidal streams, one can most re- 
liably know their positions and velocities in the past. This 
is because stars in tidal streams all originate, at one time 
in the past, from a single progenitor on a specific orbit. 
It is this property of tidal streams th at makes them excel- 
lent tracers of the Galact ic potential ( [Johnston et al.lll999l : 
iMuraU fc Dubinskil [l999l ) while the degree of tidal-stream 
coherence provides a po werful constraint on th e lumpiness 
of the Galactic potential (jjohnston et al.ll2002h . 

The disruption timescales of merging galaxies 
l|Helmi fc Whitel 1 19991 . hereafter H99a) means the dis- 
covery of new t idal s tream s, e.g. the Orph an stream 
l|Belokurov et all |2006| . l2007l : iGrillmaiij [ 20061 ). can also 
lead to the discovery of new sa t ellites, e.g. Urs a Majo r 

II (UMa II, IZucker et all l2006l : iFellhauer et"al] l2007al ). 
which are of relevance to the ' missing satellite problem ' 
l|Moore et al.lll999l : iKlvpin et al.|[l999l : iBenson et al.ll2002l ). 
These discoveries offer more evidence of the hierarchical 
structure formation of the Milky Way Galaxy and provide 
accr etion history constraints for simulations of this process 
(e.g. lAbadi et al.][200a l . If a tidal stream containing DM is 
passing through the solar neighbourhood, it would provide 
a 'cold' flow of DM particles through the numerous direct 
detection experiments on Ea rth, increasing the possibility 
of a positive direct detection l|Freese et al.ll20o3 ). 

The first evidence that there may be tidal streams in 
the solar neighbou rhood was the discovery of two streams 
l|Helmi et al.l 19991 . hereaft er H99b) in kinematic surveys 
of the local stellar halo (|Beers fc Sommer- LarsenI 1 19951 : 
IChiba fc Yoshiilll998tl. Th i s disc overy was subseque i itly co n- 
firmed bv lChiba fc Beerj l|200Ci ) and lKeplev et al.l l|2007l ). 

The most spectacular example of tidal streams is from 
the disruption of the Sagittarius (Sgr ) dwarf spheroida l 
(dSph) galaxy, originally discovered bv llbata et al.l l|l994l ). 
Estimates of the central mass-to-light ratios of dwarf 
galaxies are high, suggesting they and plausibly their 
tidal streams are all DM-dominated. The most extensive 
sky panorama of the Sgr stream is traced by M giants 
l|Maiewski et al.l 120031 ') from the Two-Micron All-Sky Sur- 
vey (2MASS). The stream is seen very clearly in the South- 
ern Galactic hemisphere, as well as heading towards the 
North Galactic Pole (NGPV iMaiewski et alT(|2003l ) found 
the stream lies along a well-defined orbital plane about 
the Galactic centre and speculated that the foreshortening 
of the stream towards the NGP means it heads back to- 
wards the Galactic plane. The Sun lies within ~1 kpc of 
tha t plane, which is within the width of the stream, lead- 
ing IMaiewski et al.l l|2003l ) to propose the possibility that 
the stream cros ses the Gala c tic pl ane near or in the solar 
neighbourhood. iFreese et al.l (|2004l ) postulated that one of 
the H99b streams is part of the Sgr stream passing through 
the solar neighbourhood and that its DM density should be 
measurable in direct detection experiments. 

iLaw et al.l (|2005l ) modelled the Sgr stream and found 



that its two arms pass through a similar position in 
the plane of the sky, near th e Virgo overdensity (VOD). 
iMartfnez-D elgado et al.l (120071 ) suggest that the VOD is part 
of the Sgr stream, which crosses the Milk y Way plane in the 
solar neighbourhood, and show that the iLaw et al] (|2005l ) 
model passes through the observed location of the VOD 
and the VOD stellar density is similar to the model pre- 
dictions. However, there is a growing weight of evidence 
against the Sgr stream existing in the solar neighbour- 
hood. Dynamical analyses of the multiple wraps of the Sgr 
stream and the newly discovered Orphan str eam found in 
the 'Field of Streams ' jSclok urov et al.ll2006l ) suggest that 



neither the Sgr 
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Newberg et al. 



Fellhauer et 



al.l l2006h 



nor Orphan streams 

.2007bl ) 3ass near the Sun. More recently, 
20071 ) used Sloan Digital Sky Survey (SDSS) 



imaging and Sloan Extension for Galactic Understanding 
and Exploration (SEGUE) spectroscopic observations of F 
turnoff and blue horizontal branch (BHB) stars to show that 
the Sgr stream is neither coincident with the VOD nor passes 
through the solar neighbourhood. Instead, they find that it 
misses the Sun by 15 kpc, passing through the plane well 
out side the sola r circle. 

lJuric et al.1 (|2005l ) used SDSS photometry to identify 
the VOD - the largest clump of tidal debris so far detected 
in the outer halo. It covers ~1000 deg^ above the Galactic 
plane (5 < ^ < 15 kpc), over the solar position in the Galaxy 
(J? ~ 7 kpc) . They i nterpreted it as a tidal stream or an 
invading dwarf galaxy. lJuric et ah! (|2005| ) did not detect any 
downturn in the star counts towards lower Galactic latitudes 
(6), indicating the VOD could extend closer to the Galactic 
plane than the observations probe {b > 60°). 

With the advent of SDSS Data Rel ease 5 and the dis - 
coveries of the VOD and Orphan stream. [Puchs et al.l (|2006l ) 
re-analysed star counts from their Calar Alto Deep Imaging 
Survey (CADIS) data. In retrospect, they find that overden- 
sities in the vertical-density distribution of stars can be as- 
sociated with t he VOD in the 13-h C ADIS field and possibly 
the Monoceros (|Newberg et al.l2002l) and Orphan streams in 
the 9-h CADIS field. Simulations bv lFellhauer et al.l (|2007al ) 
show the 9-h field falls on the second wrap of the backward 
orbit of UMa II before it passes along t he Orphan stream. 
In this direction. IFellhauer et all ([20073) demonstrated that 
the Orphan stream is >30 kpc from the Sun and its orbit 
does not bring it much closer to the solar neighbourhood 
(see their fig. 2). 

The Monoceros stream is seen in the data of lJuric et al.l 

("20051) at 3 < Z < 5 kpc and 7? ~ 16 kpc. Two other over- 
dcnsities are also visible in the top right plot of lJuric et al.l 
(2005) fig. 20 at {R, Z) ~ (1.5, 6.5) and (9, 0.8) kpc (the 
latter is seen most clearly in the second row, right plot). 
None of these four overdensities are seen to exist in the so- 
lar neighbourhood. iPefiarrubia et al.1 (|2005l ) have modelled 
the Monoceros stream and their simulations show that part 
of it may be seen in the direction of the 9-h field at heliocen- 
tric distances of the overdensity in the 9-h field (3-21 kpc), 
again far from the solar neighbourhood. 

Unlike in the Galactic disc, tidal streams in the outer 
halo are not phase-mixed quickly due to its dynamical time- 
scale of >1 Gyr and so remain coherent in configuration 
space for longer than in the disc. The low stellar density in 
the outer halo enables large tidal streams to be identified as 
stellar overdensities in photometric surveys like the SDSS. 
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The shorter dynamical timescales and much higher stellar 
densities cause tidal streams in the inner Galaxy to quickly 
lose their coherence in configuration space but they remain 
coherent in velocity space (H99a). 

Until recently, the full six-dimensional phase space of 
the solar neighbourhood had not been systematically sur- 
veyed. This was b ecause, although the Hipparcos satellite 
mission (|ESAlll997l ) provided accurate parallaxes and proper 
motions for ~118 000 stars, these stars generally lacked 
radial velocities (RVs). Therefore, a large European con- 
sortium obtained kinematically unbiased RVs of Hippar- 
cos s tars of spectral type later than about F5 (|Udrv et alj 
1 19971 ). Multi-epoch RVs of ~45 000 stars were measured with 
the two CORrelation RAdial VELocities (C ORAVEL) pho- 
toelectric cross-c orrelation spectrometers (|Baranne et al.l 
ll979l : lMavoilll98^ ). There have been two public data releases 
and three published analyses of stars in the CORAVEL 
database: the Geneva-Copenhagen survey of 16 682 nearby 
CORAVEL F-G dwarfs, available via VizieR, is described 
and analysed in Nordstrom ct al. ( 2 004) , herea fter N04, and 
recently re-analysed in IHolmberg etaU l|2007h : and a cata- 
logue of the local kinematics of 6691 CORAVEL K-M gi- 
ants, also avai l able v ia VizieR, is described and analysed by 
iFamaev etlll (120051 ). hereafter F05. 

The RAdial Velocity Experiment (RAVE, 
ISteinmetz et al.ll2006l ) is a spectroscopic survey measuring 
the RVs and stellar atmosphere parameters (temperature, 
metallicity and surface gravity) of up to one million stars 
in the range of magnitudes probing Galactic scales between 
the very local CORAVEL RV surveys and the more distant 
SEGUE RV survey. RAVE started in 2003, using the 
Six-Degree Field (6dF) multi-object spectrograph on the 
1.2-m UK Schmidt Telescope of the Anglo- Australian 
Observatory in Australia. 

The recent, timely availability of the CORAVEL and 
RAVE data sets permit direct searches for the presence o f 
tidal streams in the solar neighbourhood, [ifelmi et al.l(|2006l ) 
have already searched the N04 CORAVEL dwarf data set for 
signatures of past accretion on to the local disc. They find 
that stars with a common progenitor should show distinct 
correlations between orbital apocentre (A), pericentre (P) 
and Z angular momentum ( Jz). In APJ-space, their analysis 
reveals a statistically significant excess of stars on orbits 
of common eccentricity, analogous to the pattern expected 
for merger debris. They identify three coherent groups with 
distinct metallicity and age distributions that they assert 
correspond to the remains of disrupted satellites. 

In this paper, we (quite literally) take an orthogonal ap- 
proach by simply posing the question: is there any net verti- 
cal flow through the solar neig hbourh o od? T his is a different 
questi on to that explored by iGouldl (|2003l ) and iBeU et al.l 
(|2003): we are not attempting to constrain the amount 
of phase-space substructure. The well-established Galactic 
stellar components in the solar neighbourhood (thin and 
thick discs and inner halo) are kinematically symmetric 
about the Galactic plane. Therefore, any net vertical flow 
could reasonably be associated with a tidal stream of stars 
with vertically coherent kinematics. Hence, within the sam- 
ple volumes of the CORAVEL and RAVE surveys, we can 
directly test whether the Sgr stream and VOD pass near the 
Sun and consequently whether direct DM detection experi- 



ents should expect enhanced modulated signals from these 
or any other streams. 

We organize the paper by analysing the data sets in or- 
der of increasing sampled volume. Ergo, we begin with the 
CORAVEL dwarfs in Section [H In Section [Ol we use this 
sample to show that the effect of binarity is negligible in 
the analysis of velocity distribution functions, which is im- 
portant for our purpose here as well as many future science 
applications of RAVE and SEGUE data. We introduce our 
statistical test to identify streams in Section [2.41 This test 
is applied to the CORAVEL giants in Section [3] and all the 
RAVE stars in Section |4] 

This is not th e first science pap er to make use of 
the RAVE data set. ISmith et al] ((2003) used high- velocity 
RAVE stars to constrain the local Galactic escape speed. 
Veltz et al. (2007) used 2MASS to select RAVE stars with 
0.5 J — Ks ^ 0.7, which corresponds to K3-K7 dwarfs and 
G3-K1 (red clump) giants, to identify kinematic and density 
discontinuities between the thin disc, thick disc and inner 
halo towards the South Galactic Pole (SGP). We introduce 
the RAVE stellar populations in Section |4T] We discuss our 
results in Section [5] before ending on an historical aside in 
Appendix A. 



2 CORAVEL DWARFS 

2.1 Derivation of the space velocities of the 
CORAVEL dwarfs 

The Geneva-Copenhagen survey is complete, all-sky, 
magnitude-limited and kinematically unbiased. Its observa- 
tional input catalogue was selected from a compilation of 
catalogues available in the literature with Stromgren uvbyf3 
photom etry of nearby F and G sta rs, mainly from the sur- 
veys bv lOlsenI l| 19831 . 1 19931 . Il994al lbl). The input catalogue 
was observed using both CORAVELs. Their fixed, late-type 
cross-correlation template spectra match the spectra of the 
majority of the input catalogue stars. The multi-epoch RVs 
(generally two or more) have a modal mean error of the 
mean RVs of 0.25 km s"^ (see N04 fig. 3). 

The vast majority of these stars have proper motions 
in the Tycho-2 catalogue l|H0g et al.ll2OO0l ). This catalogue 
was constructed by combining the Tycho star-mapper mea- 
surements of the Hipparcos satellite with the Astrographic 
Catalogue based on measurements in the Carte du Ciel and 
other ground-based catalogues. The typical mean error in 
the total proper motion vector is 1.8 mas yr~^. 

The primary source of distance for these stars is Hip- 
parcos trigonometric parallax. This is adopted if its relative 
error (ct^/tt) is accurate to 13 per cent or better, otherwise 
the photometric distance calibrations for F and G dwarfs 
bvlCrawford (1975) and Olson (19M) 

are used, with an un- 
certainty of only 13 per cent. Distances are not provided 
for stars with unreliable Hipparcos parallax {a-n/n > 13 per 
cent) and when photometric distances cannot be calculated. 
This occurs when the star is missing the necessary photome- 
try and/or it falls outside the photometric calibrations. The 
absence of a distance estimate or RV measurement reduces 
the size of the sample with all full six dimensional phase- 
space information to 13 240 stars. The space velocity com- 
ponents in the Galactic cardinal directions, U (towards the 
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Figure 1. Distribution of CORAVEL dwarfs in tlie plane 
of orbital angular momentum components, where Jpiane = 

Jjr + Jy ■ Top: The triangle is star HD 175305 and the square 
is CD —80 328. Bottom: The red dots arc stars in an apparent 
linear structure. 



Galactic centre), V (in the direction of Galactic rotation) 
and W (towards the NGP), are computed for all the stars 
with (mean) RVs, proper motions and distances. 



2.2 Orbital angular momenta of the CORAVEL 
dwarfs 

The availability of distances to each CORAVEL dwarf allows 
their Galactic-centred cartesian co-ordinates {X, Y, Z) to 
be calculated, where {Xq, Yq, Zq) = (8, 0, 0) kpc, which in 
combination with their (U, V, W) relative to the Galactic 
Standard of Rest (GSR) permits their components of orbital 
angular momentum (J) to be resolved into 



Jx = YWgsr - ZVasR, 
Jy = ZUgsr — XWasR, 
Jz = XVasR — YUasR, 

where 

Ugsr = Ulsr = U + UfsR, 

VasR = Vlsr + Vrot = V + Vfsfl. + Kot, 

Wgsr = Wlsr + WfsR, 



(1) 
(2) 
(3) 

(4) 
(5) 
(6) 



where (U?,c.„, 

(IDehnen fc BinnevI 



T/G 



W 



199i 



sr) = (10.00, 5.25, 7.17) km s"' 
is the solar motion, decomposed 
into its cardinal directions, relative to the Local Standard 
of Rest (LSR), and Vrot = 220 km s~^ is the amplitude of 
Galactic rotati on towards I = 90° and h = 0° (lAU 1985 
convention; see iKerr fc Lvnden-Bell|[l986l '). Fig. [T] shows the 
CORAVEL dwarf sample is dominated by high-angular mo- 
mentum disc stars with only ~100 low-angular momentum 
halo stars. 

Two stars on prograde high-inclination orbits (the tri- 
angle and square in the top plot of Fig. [T} share simi- 
lar angular momenta as the H99b streams (cf. their fig. 
2). The triangle in the top plot of Fig. [1] is HD 175305 
and has a photometric mct allic i ty of —1.3 9 dex. It is also 
m the|B eers fc Sommer-Larsenl l|l995h and IChiba fc Yoshiil 
(I1998h catalogues with spectroscopic metallicities of —1.42 
and —1.54 dex respectively. Its kinematics, {U, V, W)gsr = 
(—58, 142, —286) km s~^, means it is on e of the stars in the 
H9 9b stream moving to wards the SGP. iBeers et all (|200G| ) 
and lKeplev et all |2003) explicitly cite this star as a member 
of one of the H99b streams. 

The square in the top plot of Fig.[T]is CD -80 328 and 
has a photometric metal li city o f —1. 98 dex. It is not in the 
iBeers fc Sommer-Larse"nl (|l995h nor IChiba fc Yoshiil (|l998l ) 
catalogue s. Its very low met allicity is spectroscopically con- 
firmed bv lBeers et al.l (|l999l ) to be —2.09 dex. Its kinemat- 
ics, {U,V,W)gsr = (-193,117,303) km s~\ and metallic- 
ity suggest it could be a member of the H99b stream moving 
towards the NGP. However, its kinetic energy is too large 
to be consistent with the ener gies of the other mem bers of 
the stream (listed in table 7 of IChiba fc Beer^l2000l and th e 
energies were calculated using table 4 of lBeers et al.ll2000l ). 
Therefore, we consider this star to be a possible outlier, like 
the different outlier found by H99b. The red dots in the bot- 
tom plot of Fig.[T]apparently align in a linear structure. The 
mode of their photometric metallicities (—0.6 dex) agrees 
with the second group of tidal debris found bv lHelmi et al.l 
l|2006t ). However, the younger and wider age distribution of 
the linear structure argues against a tidal origin, suggesting 
a chance alignment is more probable. 

The lack of strong clumping and apparent smoothness 
of the distribution in Fig. [T] suggests there are not any sig- 
nificant, coherent tidal streams in t he sample (the tidal de- 
bris identified bv lHelmi et al.ll2006l does not manifes t itself 
in coherent tidal streams) . The 'Helmi et al] (|2006l ) APJ- 
space analysis of the CORAVEL dwarfs is similar to the 
H99a entropy technique (also used in H99b), since it parti- 
tions the parameter space and counts the number of stars 
in each cell. Therefore, we do not re-apply these techniques 
to the CORAVEL dwarfs in Fig. |T] In addition, these tech- 
niques cannot be applied to RAVE data as so far we have 
not derived accurate distances to RAVE stars (see Section 
13]). Instead, we apply the RAVE data analysis technique em- 
ployed in Section |4] to the CORAVEL samples, in order to 
characterize the solar neighbourhood W distribution and to 
allow comparison between the CORAVEL and RAVE re- 
sults. First, the next section demonstrates that RAVE data 
can be used on an equal footing with the CORAVEL data. 
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Figure 2. Top: W distribution of the CORAVEL dwarfs: all the 
13 240 stars (solid histogram) and the 8589 single stars (dotted 
histogram). The 3 km velocity bin sizes are chosen to be twice 
the average space velocity error. Bottom: Normalized CDF of all 
the 13 240 CORAVEL giants (solid line) and the 8589 single stars 
(dotted line) as a function of W , where the maximum differences 
between them are indicated by vertical lines: (dot-dashed 
line) and D_ (dashed line). 



2.3 Effect of Binarity on W velocity distribution 

A double correlation peak may identify a spectroscopic bi- 
nary from a single RV measurement but an average of more 
than four Geneva-Copenhagen observations are available per 
CORAVEL dwarf, aUowing N04 to identify 3223 stars (out 
of 16 682, 19 per cent) as spectroscopic binaries of all kinds. 
If the velocities can be properly assigned to the two bi- 
nary components, the centre-of-mass velocities of double- 
lin ed spectroscop ic binaries can be computed by the method 
of IWilsonI (|l94ll ) without a full orbital solution. N04 pro- 
vide systemic velocities of double-lined spectroscopic bina- 
ries where possible, otherwise the raw average RV is given 
for CORAVEL dwarfs, including single-lined spectroscopic 
binaries. 

The average scatter of individual RVs of single-lined 
spectroscopic binaries, which have orbital periods of years 
or less, is ~15 km s~^. This scatter is an order of magnitude 
larger than the average error in each space velocity compo- 
nent of 1.5 km . Therefore, the mean RV of a single-lined 
spectroscopic binary may be more representative of the in- 



stantaneous radial component of the orbital motion in the 
binary system than its systemic line-of-sight motion in the 
Galaxy, resulting in inaccurate and even potentially mislead- 
ing space velocities. 

Although any individual space velocity could be unre- 
liable due to binarity, the statistical effect of binarity on 
the velocity distribution function of a large number of stars 
has not been empirically constrained before. This is due to 
the sparsity of RV data available with sufficient temporal 
coverage to generate good binarity detection statistics. The 
Geneva-Copenhagen survey is the largest, homogeneous RV 
survey with repeat measurements: 62 993 new RV observa- 
tions of 13 464 programme stars, representing ~1000 nights 
of data. Observations typically cover a time span of 1-3 years 
but some extend over a decade or more, allowing the major- 
ity of binaries to be identified. 

The statistical effect of binarity on the velocity distri- 
bution function is investigated for the first time here be- 
cause it has fundamental implications for our aims in this 
paper, as well as the observational strategies and science 
of on-going RV surveys. The vast majority of RAVE (and 
SEGUE) stars have never been spectroscopically observed 
before and RAVE (and SEGUE) observes most of them only 
once. Ideally kinematic studies should use centre-of-mass ve- 
locities but when these are not available, it is important to 
know how sensitive the analysis is to the kinematic affects 
of binarity. 

We investigate the kinematic effect of binarity by 
comparing the W velocity distribution function of all the 
CORAVEL dwarfs, including binaries, to the corresponding 
distribution function of single stars (see top of Fig. [2]). N04 
provide a catalogue flag (/;,) that identifies confirmed or sus- 
pected binaries, where the information can come from one or 
several sources such as photometry, RV or astrometry. They 
identified 3537 visual binaries (out of 16 682 stars, 21 per 
cent) and the majority of systems with periods below 1000 
days should be flagged. Single stars, for which the derived 
space velocities are more reliable, are defined by N04 as hav- 
ing a null /(, catalogue entry (8589 stars out of the 13 240 
with space velocities, 65 per cent). Some CORAVEL dwarfs 
only have two observations. Accordingly, some long-period 
and /or low-amplitude binary stars will not be identified by 
the fb flag and will thus be present in the sample of single 
stars. Their effect should be negligible compared to a bona 
fide sample of single stars due to the small amplitude of their 
orbi tal velocities. 

iKepley et al. I l|2007l ) were able to use a parametric test 
(the Shapiro- Wilk test for deviations from normality) to 
search for streams in their stellar halo sample because the 
halo velocity distribution is approximated by a Gaussian. 
Both W velocity distribution functions in the top of Fig. [2] 
approximately consist of the sum of three Gaussians. These 
represent the young thin disc, where aw ~ 10-15 km s~^, the 
old thin disc, where aw ~ 15-20 km s~^, and the thic k disc, 
where aw > 30 km s'^ l|Seabroke fc Gilmord [20071 ) . The 
shape of the resulting distribution is a Gaussian with posi- 
tive kurtosis (leptokurtic) : it is more peaked and has heavier 
tails than a single Gaussian. This means a non-parametric 
statistical comparison test is required. 

A standard distribution-free test is the Kolmogorov- 
Smirov (K-S) test. The K-S statistic is the maximum dif- 
ference over all values of a single, independent variable x 
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of two cumulative distribution functions (CDFs). However, 
the sensitivity of the K-S test is not independent of a;: it 
tends to be most sensitive around the median value and less 
sensitive at the extreme ends of the distribution. Identify- 
ing the points at ±00 (wrapping the a;-axis around a circle) 
guarantees equal sensitivities at all valu es of x. This is a n 
invariant K-S test called the Kuiper test (|Press et al.ll 1993 ). 
For comparing two different CDFs, the Kuiper statistic, de- 
fined as 

D = D+ + D., (7) 
is the sum of the differences between the CDFs, where 
D+ = max[5'jvi (x) — S'jv2(a;)] (8) 
and 



D. 



max[SN2 (^) ~ 'S'iVi (x)], 



(9) 



where SNi{x) is the function giving the fraction of data 
points ^x. The statistical significance of D, P{D > 
observed) = Q{\) is given by 



Q = 2^(4/A^-l)e-^^"^^ 



(10) 



where Q is a monotonic formula with variable j for the 
asymptotic behaviour of the statistic D, satisfying Q{0) = 1 
and Q{oo) — 0, where 



X = D 
and 



/A^e +0.155 + 



0.24 

7^ 



iVl + iV2 ' 



(11) 



(12) 



To compare all the CORAVEL dwarfs (iVi = 13 240) 
with the single-star CORAVEL dwarfs {N2 = 8589), x = 
W in equations M and H giving D+ = 0.004, D- = 0.009 
(see bottom of Fig. [2]) and D = 0.013. The sample size is 
sufficiently large (A^e = 5209.51) to give the test enough 
power to detect significant differences and reject the null 
hypothesis (that all the stars and the single stars are drawn 
from the same parent population) if it is false. Q = 0.887, 
signifying the data fail to reject the null hypothesis. This 
means that, statistically, stellar binarity does not affect the 
W distribution of the CORAVEL dwarfs significantly. 

Because binaries are, on average, apparently brighter 
than single stars, binaries are abundant in apparent magni- 
tude limited samples like the Geneva-Copenhagen survey. 
The Kuiper test suggests we can safely include binaries, 
nearly doubling the sample size of the CORAVEL dwarfs 
and increasing the power of the test. This means the RAVE 
RVs, which includes binaries, can be used on an equal foot- 
ing with the CORAVEL data. Fig. [3] shows that including 
binaries does not increase the local volume of the solar neigh- 
bourhood sampled by the CORAVEL dwarfs. 



2.4 Determining Wlsr for the CORAVEL dwarfs 

To search for streams, a Kuiper test could be used to com- 
pare the W distribution of the CORAVEL dwarfs with 
a Cauchy-Lorentz distribution (heavy-tailed Gaussian). Its 
peak location (a;o) and half-width at half-maximum scale 




Figure 3. Three dimensional sky distribution of all the 13 240 
CORAVEL dwarfs with W velocities (red) and the 8589 single- 
star CORAVEL dwarfs (blue). The spherical polar axes have a 
radius of 370 pc, centred on the Sun, chosen to include the most 
distant CORAVEL dwarf with a W velocity (363 pc, marked by 
the black cursor). The viewer is in the Galactic plane, inside the 
solar circle at I ~ 30° (the Galactic centre is over the viewer's 
left shoulder). The nearest line of constant I is I = 30° (with the 
majority of the stars behind it). I increases to the right (anticlock- 
wise) with lines of constant I every 30° . All the lines of constant I 
converge at the NGP (b = +90°, top) and at the SGP {b = -90°, 
bottom) , with lines of constant b every 30° . 



(7) parameters could be theoretically determined accord- 
ing to the Galactic volume sampled and the age-velocity 
dispersion relation of the thin and thick discs. However, 
the functional form of th ese relations remains uncertain 
(Sea broke fc ^ilmoro"2007|). In addition, the test hypothesis 
would be whether the data is exactly a Cauchy-Lorentz dis- 
tribution with the specified xo and 7. The hypothesis could 
be rejected due to the data having a slightly different xo 
and 7 than those specified. Thus, the test may not be very 
sensitive to the presence of streams. 

A more straight-forward approach is to test the symme- 
try between the two halves of the W distribution. In-falling 
streams on to the Milky Way disc will reveal themselves 
as overdensities in the distribution of W with respect to the 
Galactic plane. W is positive towards the NGP and negative 
towards the SGP (symbolically denoted +W and —W re- 
spectively for brevity), so a stream falling onto the Galactic 
disc from the NGP ('above' the Galactic plane) would have 
~W and a stream falling onto the Galactic disc from the 
SGP ('below' the plane) would have +W. A stream falling 
through the disc from the NGP to SGP would exhibit ~W 
both above and below the plane and +W both above and 
below the plane falling from the SGP to the NGP. A sym- 
metry test between the two sides of the W distribution, if 
sensitive enough to the number of stream stars, will find the 
asymmetry caused by a single stream. 

W is measured relative to the Sun. The top of Fig. [2] 
shows that the W distribution is not centred on zero velocity. 
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IWlsrI (km s-') 

Figure 4. Top: Wlsr distribution of the CORAVEL dwarfs 
(solid histogram) either side of the Wlsr = dividing line 
(dashed line). Bottom: Normalized CDF of W^sR < (solid line) 
and Wlsr > (dotted line), showing D— (dashed line) and 
(dot-dashed line). 

The Kuiper test is sensitive to the sign of the W velocities 
(see bottom of Fig. \^ so the absolute values oi — W veloc- 
ities are required to compare them to -\-W to test the sym- 
metry of the distribution. The easiest method to implement 
this Kuiper symmetry test is to transform the heliocentric 
reference frame to one centred on zero velocity - the LSR 
reference frame, where in equation |6] is the centre of 

the W distribution in the top of Fig. [2] W^fg^ needs to be 
determined to apply its correction to the W distribution to 
convert it to a Wlsr distribution. 

We determine W^^^ by varying its value in the Kuiper 
symmetry test to find the best agreement between the 
two halves of the resulting Wlsr distribution (minimum 
value of D) . D is calculated using equations [8] and |9] but 
now X = |±M^LSfl|- The maximum differences between 
'S'iVw^g^<o(|-WLSfll) and S'jv„.^^,^^o (+Wlsh) are calcu- 
lated over the range < |±VyLSii| < oo. The top of Fig. 
|4] shows the Wlsr distribution for the minimum value of 
D, where W^g^ = 7.0 km s~^ to the nearest 0.1 km s~^. 
Our value of W^ap, also agrees within the errors with the 
iDehnen fc Binnevl lfl993) measurement: W^g^ = 7.17 ± 
0.38 km s~^ (see Appendix A for further discussion on the 
determination of W®g^ and historical parallels). 

Having determined Wfg^, we continue to vary its value 



until a 4a" statistically significant difference (Q < 6 x 10~^) 
is generated between the two halves of the resulting Wlsr 
distribution to gauge how sensitive the Kuiper symmetry 
test is to the measured value of W^g^^. Table 1 shows that 
'^LSR '^^i placed —1.4 km s~^ from our measured value 
of in the —Wlsr distribution, whereas it can only 

be placed less than half that speed (-1-0.6 km s~^) from 
our measured value of Wfg^ into the +Wlsr distribution. 
This is because the position of W^g^ is sensitive to the 
slight asymmetry due to D+ in the bottom of Fig. |4] at 
IWLSfil ~ 13 km s~^. Although it may not appear so in the 
top of Fig. m because of the binning, there are more stars 
with < Wlsr sS 13 than -13 < Wlsr < km s^\ How- 
ever, Sn^.^^^^oH-WlsrI) > SiVK^^^^^J+WiSii), allowing 
the line of symmetry (W^g^) to be placed further into the 
— Wlsr distribution than the +Wlsr distribution. If the 
position of W^gj^ was sensitive to the larger asymmetry due 
to D- in the bottom of Fig. g] at |Wz,sfl| = 25 km s"\ 
the reverse would be true i.e. even though there are more 
stars with < Wlsr ^ 25 than -25 < Wlsr sS km s"\ 
SN„^^^>oi+'^LSR) > Snw^^^^oH-Wlsr]), which would 
allow the line of symmetry (W^g^) to be placed further into 
the +Wlsr distribution than the —Wlsr distribution. This 
demonstrates that is only sensitive to asymmetries at 

IWLSfil < f^w^sR (in i-his case -D+ rather than D-), despite 
D+ « D„. 

'^LSR °^ tii6 full sample of CORAVEL dwarfs coin- 
cides with the mode of the Wlsr distribution in the top of 
Fig. |4l This suggests that there cannot be many CORAVEL 
dwarfs with < cwlsr belonging to vertically coher- 

ent streams with systemic Wlsr velocities. If there were, 
they would have biased the measurement of away 
from the mode. 

A coherent stream with a systemic space velocity would 
be visible as an overdensity in the CORAVEL dwarf U—W 
and V—W space diagrams. These diagrams in fig. 20 in 
N04 do not exh ibit any such vertical substructure. In fact, 
ISeabroke fc~Gil morc (2007) show the CORAVEL dwarfs are 
well phase-mixed vertically because of the sample's low 
scaleheight (see Fig. [Sjl . This suggests there are no vertically 
coherent streams currently passing through the solar neigh- 
bourhood, in agreement with the lack of strong clumping in 

Fig.m 



2.5 Sensitivity of the Wlsr symmetry of the 
CORAVEL dwarfs to tidal streams 

In this section, we investigate the sensitivity of the Wlsr 
symmetry of the CORAVEL dwarfs to tidal streams at 
VKlshI ^ ""Wj^sH- The previous section suggests the sam- 
ple of CORAVEL dwarfs is free from vertically coherent 
streams. This means its distribution can be used as a smooth 
background against which we can investigate the effect of 
a tidal stream passing vertically through the solar neigh- 
bourhood. We do this by simulating a stream, pseudo- 
randomly generating a Gaussian with as = 10 km s~^, 
and adding it to the sample. This dispersion was chosen 
because the best studied tidal stream in the literature, the 
Sgr stream, is found to be dynam ically cold by two differ- 
ent groups: iMaiewski et al.l l|2004l ) used medium-resolution 
spectroscopy to find an intrinsic dispersion of a = 10.4 ± 1.3 
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Table 1. Results of the Kuiper symmetry test, applied to the CORAVEL dwarfs, to determine Wj^gj^ to the nearest 0.1 km s ^, 
where AVF®g^ = VKf^^ - 7.0 km s"!. 



WisR 






Nwlsr>0 




D_ 


D+ 


D 


Q 


5.6 


-1.4 


7176 


6064 


3287 


0.0484 


0.0021 


0.0505 


3 X 10-*^ 


7.0 


0.0 


6583 


6657 


3310 


0.0250 


0.0064 


0.0314 


0.0341 


7.6 


0.6 


6334 


6906 


3304 


0.0242 


0.0212 


0.0454 


6 X 10-5 



km s ^ (after r emov ing random errors of ~5.3 km s ^); and 
iMonaco etHI (|2007l ') fitted a Gaussian of a = 8.3 ± 0.9 km 
s"^ to their sample of high-resolution spectroscopically de- 
rived RVs. 

To investigate how sensitive the Kuiper symmetry test 
is to the choice of as , we repeated the above experiment with 
o"s = 3 and 17.3 km s~^. These dispersions are physically 
motivated: H99a state that the intrinsic dispersion for a tidal 
stream from a Large Magellanic Clou d-type progenitor is 
~3-5 km s"^ after a Hubble time; and lDuffau et af] (|2006l ) 
discovered the Virgo stellar stream with a = 17.3 km . 
We find that the number of stars required to generate a 
4a tidal stream detection is relatively insensitive to as- The 
conservation of phase-space density results in a decreasing 
as 1/time at each point along a stream (H99a). Hence, the 
ability of the test to detect a stream is not biased by the 
dynamical age of the stream. 

The simulated tidal stream is added to the sample of 
CORAVEL dwarfs at W = ±(1, 2, 3, ^13)awLSR' where 
awi^sR ~ ^8 t^^* sensitivity of the Wlsr 

symmetry to the presence of each stream. Approximately 
— ISa wr va corresponds to Wlsr ~ —242 km calcu- 
lated by Freese et al. ( 20041 ) using the eight clump stars from 
IChiba fc Yoshiil (|l998( ) i"nthe H99b stream (the CORAVEL 
dwarfs with Wlsr nearest to this value are at —286 and 
-179 km s-i). 

At each position, the random population of the sim- 
ulated tidal stream is increased until it causes the Kuiper 
symmetry test to reject the null hypothesis (that Wlsr < 
and Wlsr > are drawn from the same parent population) 
at the 4cr statistical significance level. The test results in 
Table 2 are indicative because they represent a single real- 
ization of a pseudo-random Gaussian added to the sample 
and we have not conducted a full Monte Carlo simulation. 
Nevertheless, the indication that more tidal stream stars are 
required to generate a 4a tidal stream detection at 2 and 
Sawj^sR than at —2 and ^3awLSR. seems plausible due to 
the asymmetry already in the Wlsr distribution: D- in the 
bottom of Fig.|4]at | WlshI = 25 km . Placing streams in 
the —Wlsr distribution, further increases D- and decreases 
D+. The reverse is true for placing streams in the +Wlsr 
distribution but even more stream stars are required to can- 
cel out the intrinsic D_ and increase D+ to the 4a statistical 
significance level (see Fig. [5} . 

The test becomes more sensitive (less stream stars re- 
quired to generate a 4a detection), the further a stream is 
from the centre of the distribution. However, Table 2 shows 
that the number of stream stars required for a 4a detection 
stays approximately constant (saturates) at ^—Sawj^sR- 
This is because the positions of D- and D+ have become 
fixed and their values saturated, whereas D+ is still sensitive 
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Figure 5. Top: Wlsr distribution of the CORAVEL dwarfs 
(solid histogram) either side of the Wlsr = dividing line 
(dashed line). The pseudo-randomly generated Gaussian (rep- 
resenting a cr = 10 km s" ^ tidal stream) included in the data 
set (dotted histogram) that caused the Kuiper symmetry test to 
reject the null hypothesis at 4cr when placed at ■^crWi^sR ^^^^ 
plotted separately (dot-dashed line). Bottom: Normalized CDF of 
Wlsr < (solid line) and Wlsr > (dotted line), showing _D_|_ 
(dot-dashed line). D— is not visible in this plot at IWlshI ~ 0. 



to perturbation at Wlsr > 3crw£,sj, by hundreds of stream 
stars. 

H99a demonstrate that phase-mixing allows different 
tidal streams from a single disrupted object to be observed 
with opposite motion in the vertical direction (e.g. the two 
H99b streams). If two such streams contain similar numbers 
of stars, the Kuiper symmetry test may not find them be- 
cause their presence may not create enough of an asymmetry 
between the two sides of the W distribution. H99a predict 
inside a 1 kpc^ volume centred on the Sun, there should be 
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Table 2. Results of the Kuiper symmetry test, applied to the CORAVEL dwarfs, to determine the approximate number of tidal 
stream stars (Ns) required, in pseudo-randomly generated as = 10 km s~^ Gaussians, placed at ±(1, 2, 3, ~13)(Tvi'j^gjj in the 
Wi^SR distribution, to cause the test to reject the null hypothesis at 4cr (Q < 6 X 10~^). % = Ns/{Ns + Ntotal) where Nt^tal = 13 
240. Ps is the stellar density of each stream, in Af X 10® stars kpc~'^, calculated by dividing the Ng of each stream by the volume 
within which the CORAVEL dwarfs are complete, approximated by a sphere with radius ~40 pc {~3xl0^* kpc"'^). 





Ns 


% 


Ps 


Nwlsr<0 


Nw^SR>0 






D+ 


D 


Q 




-1 


600 


4.3 


2.2 


7165 


6675 


3456 


0.0414 


0.0033 


0.0447 


5 X 10" 


-6 


1 


500 


3.6 


1.9 


6601 


7139 


3430 


0.0163 


0.0339 


0.0502 


2 X 10" 


-6 


-2 


200 


1.5 


0.7 


6783 


6657 


3360 


0.0469 


0.0010 


0.0479 


1 X 10^ 


-5 


2 


400 


2.9 


1.5 


6583 


7057 


3406 


0.0026 


0.0458 


0.0484 


7 X 10" 


-6 


-3 


200 


1.5 


0.7 


6783 


6657 


3360 


0.0508 


0.0004 


0.0512 


1 X 10" 


-6 


3 


400 


2.9 


1.5 


6583 


7057 


3406 


0.0026 


0.0504 


0.0530 


3 X 10" 


-7 


~-13 


200 


1.5 


0.7 


6783 


6657 


3360 


0.0508 


0.0004 


0.0512 


1 X 10" 


-6 


-13 


300 


2.2 


1.1 


6583 


6957 


3382 


0.0036 


0.0445 


0.0481 


9 X 10" 


-6 



300-500 tidal streams containing a few hundred stars each if 
the whole stellar halo was built up by disrupting satellites. 
However, neither the CORAVEL nor RAVE surveys provide 
complete coverage of this volume and so are very incomplete 
with respect to halo stars. Thus, if these streams exist in our 
samples, they will most likely contain fewer stars than the 
detection thresholds of the Kuiper symmetry test in Table 
2. This means the test is not sensitive to these streams, re- 
gardless of their possible symmetries. The number of stars 
required for a detection suggests the test is better suited to 
detect less phase-mixed, more coherent streams in configu- 
ration space, like the Sgr streams that contain hundreds of 
stars. 

'Seabroke fc Gilmord (|20o3) estimated that are ^1300 
CORAVEL dwarfs in the Hercules dynamical stream. This 
substructure corresponds to 10 per cent of the sample and 
is easily visible in the N04 U—V space diagram (their fig. 
20). F05 show in their sample of CORAVEL giants that the 
Hercules stream has au ~ 28.35±1.68 and ay ~ 9.31 ±1.22 
km s~^. More than halving au would also more than halve 
the number of Hercules stars, bringing them into agreement 
with the numbers in Table 2 but its nucleus would still be 
discernible as an overdensity in the N04 U—V space dia- 
gram. 

The high number of stream stars needed to cause a sta- 
tistically significant asymmetry would most likely be visible 
as substructure in the N04 U—W and V—W space diagrams 
and strong clumping in Fig. [T] Given this and the symme- 
try of the Wlsr distribution, we conclude that it is unlikely 
that the sample of CORAVEL dwarfs contains any verti- 
cally coherent streams consisting of >200 stars. Veltz et al. 
(2007) find the scaleheight of the thin disc is ~225 pc so the 
CORAVEL dwarfs only sample a very local volume of the 
thin disc (N04 showed the local thick disc consists of <10 
per cent of the sample). 

F-G turnoff stars are a good indicator of stellar number 
density. The only part of the Sgr stream close enough to the 
Sun where this quantity can be relia bly measured is 50 kpc 
from the Sgr dSph (SDSS Stripe 82'). lFreese erall (|2005l ) es- 
timate the full- width at half-maximum of the Sgr stream in 
Stripe 82 is 4-8 kpc, which is consistent wit h o" ~ 2 kpc, de- 
rived from 2MASS M gia nts in the stream l|Maiewski et al.l 
l2003l ). lFreese et al.l l|2005^ use this stream width to estimate 
its stellar density: 210-740 stars kpc~^. They use Blue Hor- 
izontal Branch stars in Stripe 82 and in another part of the 



stream, also 50 kpc from the Sgr dSph (Stripe 10), to es- 
timate how stellar density varies as a function of position 
along the stream and find the variation is a factor of 1-2, 
implying the stellar density of the Sgr stream in the solar 
neighbourhood is in the range of 210-1480 stars kpc"^. An 
upper limit on the stream density of 1.6 x lO'^ Mq kpc~^ is 
derive d from the models shown in iMarti'nez-Delgado et al.l 
(|2007l ) , assuming all the mass in the Sgr dSph is in the form 
of stars. 

The CORAVEL dwarfs are volume complete out to ~40 
pc, which corresponds to a spherical volume of ~3 xlO"'* 
kpc'^. Therefore, even the upper bound on the number of Sgr 
stream stars expected to be passing through this volume is 
<1. Fig. 13] illustrates that there are a significant number 
of stars out to ~ 100-200 pc so if the Sgr stream is pass- 
ing through the spherical CORAVEL dwarf volume (radius 
~200 pc), only of order unity CORAVEL dwarfs are likely 
to be Sgr stars. Table 2 shows that the Kuiper symmetry 
test is not sensitive to such a small number of stream stars 
so if there is subtle clumping in Fig. [T] the test would not 
find such weak str eams. 

From fig. 1 of lFuchs etakl l|200d ). we estimate that at 
Z ~ 2 kpc, the VOD minus the vertical density profile at 
the position of the Sun of a smooth Galaxy model is 1.2 
X 10^ stars kpc'^* in the 13-h CADIS field. Therefore, if this 
structure exists in the solar neighbourhood with this density, 
there will not be many more than 30-40 CORAVEL dwarfs 
belonging to it. Table 2 shows that the Kuiper symmetry test 
is not sensitive to such a small number of stars. The table 
illustrates that even the smallest number of stars required 
for a stream detection in this sample only puts a very weak 
constraint on the stellar density of a coherent structure in 
the solar neighbourhood with vertically coherent kinematics. 



3 CORAVEL GIANTS 

3.1 Derivation of the space velocities of the 
CORAVEL giants 

The F05 sample of CORAVEL giants is the intersection of 
several data sets: (i) spectral types K and M in the Hipparcos 
catalogue; (ii) proper motions from the Tycho-2 catalogue; 
(in) CORAVEL RVs for stars in list (i) in the Northern 
hemisphere {S > 0°), observed with the Swiss 1-m telescope 
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at Observatoire de Haute-Provence, France. The Hippar- 
cos parallaxes were used to construct a crude Hertzsprung- 
Russell (H-R) diagram from which dwarfs were identified 
and removed from the sample. Stars with peculiar spec- 
tra, such as T Tau stars, Mira variables and S stars were 
eliminated using diagnostics based on RV variability and 
CORAVEL cross-correlation profiles, combined with litera- 
ture searches. Binaries for which no centre-of-mass RV could 
be estimated were removed using the F05 binarity fiag (B). 
A star was excluded from their (and our) kinematic analysis 
if it is a spectroscopic binary {B = 0) or a visual binary 
with no orbit available {B — 5) or an uncertain case, either 
a spectroscopic binary or supergiant {B = 3) or for a vi- 
sual binary {B = 8), leaving 6030 stars (5311 K and 719 M 
giants) . 

Because the CORAVEL giants are at distances approx- 
imately an order of magnitude greater than the CORAVEL 
dwarfs, the aTr/ir distribution for the giants is correspond- 
ingly less accurate than that for the dwarfs. This means, 
unlike for the dwarfs, simple inversion of the parallaxes of 
the giants results in biased estimates of their distances for 
the > 3000 giants with aTr/ir > 20 per cent (Brown ct al. 
liggTl ). To derive bias-free distances for all the CORAVEL 
giants, F05 used t he Bayesian Lur i-Mennessier maximum- 
likelihood method l|Luri et ahll 19961 ) to model the sample in 
order to assign each star to a kinematic 'base' group. Given 
the maximum-likelihood estimator of its assigned group pa- 
rameters and its observed values, F05 obtain the marginal 
probability density law for the distance of each star from 
the global probability density function. These distances were 
used with the CORAVEL RVs and Tycho-2 proper motions 
to derive the space velocities. 

3.2 Orbital angular momenta of the CORAVEL 
giants 

Fig. [6]shows that, despite probing a larger volume than the 
CORAVEL dwarfs, the giants, unlike the dwarfs, do not in- 
clude a symmetric (no net rotation) halo sample. This figure 
includes all the sample and thus none of the CORAVEL gi- 
ants are members of the H99b streams. The FQ5 dynamical 
streams are part of the smooth distribution in Fig. [S] because 
all the stars exhibit a similar range of thin disc Jz values. 
The lack of strong clumping again suggests there are not 
any coherent tidal streams in the sample. 

3.3 Determining Wlsb^ for the CORAVEL giants 

This section repeats the technique applied to the CORAVEL 
dwarfs in Section [2]4] to the CORAVEL giants. The top of 
Fig.[7]shows the Wlsr distribution for the minimum value of 
D, where = 7.0 km s~^ to the nearest 0.1 km s~^ . This 

is the same value as found for the CORAVEL dwarfs. In the 
reverse situation to the dwarfs, VK^g^ for the giants can be 
placed further in the +Wlsr distribution than the —Wlsr 
distribution. The placement is twice as far as the maximum 
^LSR position for the dwarfs. This shows the position of 
^LSR the giants is now sensitive to the slight asymmetry 
due to _D_ in the bottom of Fig. [7] at |W^lsh| = 24 km s~^ 
because D+ is at the centre of the distribution. 

^LSR ill agreement with the mode of the Wlsr dis- 
tribution in Fig. [T] Again, this suggests that there cannot 
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Figure 6. Distribution of CORAVEL giants in the plane of or- 
bital angular momentum components, colour-coded according to 
the F05 maximum likelihood base group assignment to known 
kinematic features of the solar neighbourhood: smooth back- 
ground (black), Sirius moving group (magenta), young kinemat- 
ics (yellow), Hyades-Pleiades supercluster (red), Hercules stream 
(green) and high- velocity stars (blue). Note the different scale 
between the top plot above and in Fig. [T] 



be many stars belonging to vertically coherent streams with 
systemic Wlsr velocities. If there were, they would have bi- 
ased the measurement of away from the mode. Table 
3 shows the value of D- for the giants is greater than for the 
dwarfs and so the value of D is also larger for the giants than 
for the dwarfs. In spite of this, the Q value for the giants is 
less significant than for the dwarfs because the A^'e value for 
the giants is half the value for the dwarfs. Hence, the D- 
asymmetry does not bias the W^g^ value for the giants. 



3.4 Sensitivity of Wlsr symmetry of the 
CORAVEL giants to tidal streams 

This section repeats the technique applied to the CORAVEL 
dwarfs in Section 12.51 to the CORAVEL giants, except the 
dispersion of the giants is slightly larger: uwlsr ~ 19 

. Table 4 shows that, because D+ is at the centre of the 
distribution, many more stream stars are required to cause 
D+ to generate a 4cr stream detection at Wlsr = ^'^w^sr 
than at —law^sR - before, the test becomes more sensi- 
tive (less stream stars required to generate a 4ct detection), 
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Table 3. Results of the Kuiper symmetry test, applied to the CORAVEL giants, to determine Wf^gji to the nearest 0.1 km s ^ , 
where AVF®g^ = ly®^^ - 7.0 km s"!. 



W^SR 




Nwlsr<0 






D- 


D+ 


D 


Q 


5.4 


-1.6 


3397 


2633 


1483 


0.0609 


0.0082 


0.0691 


3 X 10-5 


7.0 


0.0 


3127 


2903 


1505 


0.0316 


0.0081 


0.0397 


0.1415 


10.4 


3.4 


2589 


3441 


1477 


0.0016 


0.0677 


0.0693 


3 X 10-5 



Table 4. Results of the Kuiper symmetry test, applied to the CORAVEL giants, to determine the approximate number of tidal 
stream stars (Ns), required in pseudo-randomly generated as = 10 km s-^ Gaussians, placed at ±(1, 2, 3, ~13)(Tw-j^g^ in the 
Wlsr distribution, to cause the test to reject the null hypothesis at 4(T (Q < 6 X 10~^). % = Ng/(Nii + N^^^^i) where N^^^^i = 
6030. ps is the stellar density of each stream, in N stars kpc— ^, calculated by dividing the Ns of each stream by the volume within 
which the CORAVEL giants are complete, approximated by a hemisphere with radius ~0.29 kpc (~0.05 kpc" '^). 
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pa 
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0.0084 


0.0688 
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16 000 
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0.0665 


0.0704 


2 X 10" 


■6 


-2 
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3327 
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0.0772 


0.0064 


0.0836 
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3127 


3203 
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0.0687 


0.0739 
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4000 


3327 
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0.0804 


0.0064 


0.0868 


5 X 10" 
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0.0797 


2 X 10" 
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200 
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4000 


3127 
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1557 
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0.0816 
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the further a stream is from the centre of the distribution. 
However, Table 4 shows that, again, the number of stream 
stars required for a Ag detection stays approximately con- 
stant (saturates) at ^— 3(Tw£,sji for the same reason as be- 
fore. 

The number of stream stars required in the dwarf sam- 
ple are also needed in the giant sample at Wlsr = ^'^Wlsr 
to generate a 4o- detection (see Fig. |8]). This is because the 
top of Fig. [5] shows that both samples include similar num- 
bers of stars at Wlsr = iSavKi^jj, even though the dwarf 
sample is twice as large as the giant one. The larger vol- 
ume sampled by the giants (see Fig. llUp increases the num- 
bers of giants in the tails of its Wlsr distribution. It is this 
greater volume, and correspondingly different stellar popu- 
lation mix, that makes the giant distribution differ from the 
dwarf one as seen in the bottom of Fig. [O] 

Fig. [TU] shows that the nine stars at ~4.55 kpc to- 
wBJds the NGP , proposed as a candidate halo stream by 
iMaiewskil (Il992l'). are not i n the sample of CORAVEL gi- 
ants. iMaiewski et al.l (|l994l ) measured the RVs of six out of 
the nine stars and H99a calculated their mean vertical ve- 
locity to be {W) = -76 ± 18 km and o-py = 35 ± 24 
km s-^. There are giants in the sample with these veloc- 
ities but they do not share the angular momenta of these 
stars computed by H99a as Jz << 0. This is far from the 
range in Fig. [S] suggesting none of the CORAVEL giants 
are members of this stream. We do not see any evidence 
of the more general halo streaming from the NGP tow ards 
the Galactic plane r e ported by[M aicwski et al.l (|l996l ) and 
iKinman et al.l l|2007l ). lKinman et al.j (i2007i ') speculated their 
streaming could be connected with the VOD. 

The method briefly described in Section 13.11 allowed 
F05 to search for and characterize kinematic substructure in 
their sample. Their method identified three kinematic base 
groups: the Sirius moving group, giants with young kinemat- 



ics and the Hyades-Pleiades supercluster. F05 found these 
groups have dispersions similar to cold streams {ctu,v,w ~ 
5-10 km s-^). The smooth background contains more than 
an order of magnitude more stars than in each group. The 
groups are barely visible as overdensit ies against the back- 
ground (see lSeabroke fc Gilmorell2007l fig. 3). Each contains 
similar numbers of stars (see table 2 in F05) to the numbers 
in Table 4 required for our Kuiper symmetry test to detect 
vertically coherent streams. F05 did not report any vertical 
substructure. This agrees with our Kuiper symmetry test 
results and t he CORAVEL g iants being well phase-mixed 
vertically (see ISeabroke fc Gilm orc 2007 figs 4 and 5). This 
suggests the D- asymmetry in Fig. [7] is not caused by a 
vertically coherent stream. 

Approximately 80 per cent of the CORAVEL giants are 
Hipparcos 'survey' stars. Therefore, the sample is complete 
for the K and M giants brighter than V = 7.3 -I- l.l|sin6| 
(|Udrvet al.lll997l ). This magnitude threshold translates into 
distances of 290 pc in the galactic plane and 480 pc in the 
direction of the galactic pole for a typical Mv = giant. 
These values are, however, very sensitive to the adopted ab- 
solute magnitude. They become 45 pc for a subgiant with 
Mv ~ +4 in the galactic plane, and 2.9 kpc for a super- 
giant with A4v ~ —5 (e.g. the most luminous giants and 
supergiants colour-coded yellow in Fig. I10|l . We consider the 
volume within which the CORAVEL giants are complete to 
be approximated by a hemisphere with radius ~0.29 kpc 
(~0.05 kpc-^). 

From the iFreese et al.l (|2005l ) arguments, we estimate 
that 10-80 Sgr stream stars could be within the volume 
where the CORAVEL giants are complete. A caveat against 
interpreting the D- asymmetry as not due to a stream is 
that the numbers of stars found by F05 to be part of sub- 
stru cture are much la rger than the expected number of Sgr 
and IMaiewskil l|l992l ) halo stream stars. This implies that 
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Figure 7. Top: Wlsr distribution of the CORAVEL giants (solid 
histogram) either side of the Wlsr = dividing Une (dashed 
hne). The 5 km velocity bin sizes are chosen to be twice 
the average space velocity error. Bottom: Normalized CDF of 
^LSR < (solid line) and Wi^sR > (dotted line), show- 
ing D— (dashed line) and D+ (dot-dashed line, just visible at 
\Wlsr\ ~ 0). 



Figure 8. Top: W^sR distribution of the CORAVEL giants (solid 
histogram) either side of the Wlsr = dividing line (dashed 
line). The pseudo-randomly generated Gaussian (representing a 
a = 10 km s~^ tidal stream) included in the data set (dotted his- 
togram) that caused the Kuiper symmetry test to reject the null 
hypothesis at 4(7 when placed at ■^owj^sr plotted sepa- 

rately (dot-dashed line). Bottom: Normalized CDF of W^sR < 
(solid line) and W^sR > (dotted line), showing D_ (dashed 
line) and (dot-dashed line). 



both the F05 method and the Kuiper symmetry test cannot 
rule out the presence of these streams in the volume sampled 
by CORAVEL giants. 

The larger volume sampled by the giants than the 
dwarfs can, however, rule out the prese nce of the VOD, i f 
their stellar densities estimated from the lFuchs et all (120061 ) 
CADIS fields are realistic. If this is the case, they would con- 
sist of 6000-8000 CORAVEL giants, which the Kuiper test 
can strongly rule out. Table 4 shows that the test can rule 
out these streams down to a density of ~4000 stars kpc""^ 
(depending on kinematics). 

Thick disc and halo stars make up <10 per cent of the 
sample (high- velocity stars, colour-coded blue in Fig. [TO)) . 
Therefore, although the CORAVEL giants have extended 
our search volume to greater vertical distances than the thin 
disc scaleheight (~220 pc), it is still dominated by the thin 
disc. RAVE is the only RV survey that vertically samples 
significant numbers of thick disc stars (as well as old thin 
disc stars) at the solar position. 



4 RAVE 

4.1 RAVE stellar populations 

To date (2007 June 26), 220 070 RAVE spectra have been 
amassed from 196 131 stars. This paper utilizes an in- 
ternal RAVE data release containing 151 856 RVs, which 
includes 24 748 RVs from the first public data releas^B 
(|Steinmetz et al.ll2006l ). a similar number of currently un- 
published RVs from the forthcoming second public data re- 
lease (Zwitter et al. 2007, in preparation) and future data 
releases (see Fig. Ill|l . 

Unlike the CORAVEL surveys that targeted specific 
stellar populations (F-G dwarfs and K-M giants), RAVE 
does not employ any colour selection or spectral type cuts. 
This means RAVE observes all the stellar populations within 



^ The first public data release of the RAVE catalogue is available 
and can be retrieved or queried from the RAVE Collaboration 
website; |http:/ /www.rave-survey.org| 
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Figure 9. Top: W^sR distribution of the CORAVEL survey: 
6030 giants (solid histogram) and 13 240 dwarfs (dotted his- 
togram). Bottom: Normalized CDF of the giants (A'^i, solid line) 
and the dwarfs (N2, dotted line) as a function of Wlsr^ where 
the maximum differences between them are indicated by vertical 
lines: D+ (= 0.0509, dot-dashed line) and D- (= 0.0257, dashed 
line), where = 4143.08, D = 0.0766 and Q = 1 X 10'^^. 

the apparent I magnitude-limit ed selection func t ions o f its 
input catalogues (see Fig. [HJ. ISteinmetz erall (|2006l ) fig. 
5 shows that this selection function is much fainter than 
the CORAVEL (and Hipparcos) surveys. Therefore, there is 
only a few per cent serendipitous overlap of RAVE stars with 
Hipparcos parallaxes (with large errors) and/or Stromgren 
photometry at the bright end of RAVE. Hence, the dis- 
tance determination methods described in Section 12.11 are 
not available to RAVE stars and the method described in 
Section 13.11 is beyond the scope of this study. RAVE stars 
cannot be placed in a H-R diagram because sufficiently accu- 
rate photometrically derived distances are not yet available 
(discussed in more detail below). 

Instead, we illustrate RAVE's different stellar popula- 
tions using a reduced proper motion diagram. The concept of 
reduced prope r motion [H) was first used by E. Hertzsprung 
l|Luvtenlll96a ) in the absence of absolute magnitude. Just 
as parallax fixes the absolute magnitude exactly, H deter- 
mines the absolute magnitude approximately. When this is 
plotted against colour, it produces a statistical H-R dia- 
gram (kinematic colour-apparent magnitude diagram). H 
combines the observable properties of proper motion {/J., in 




Figure 10. Three dimensional Galactic sky distribution of all 
the 13 240 CORAVEL dwarfs with W space velocities (grey) and 
the 6030 CORAVEL giants with W space velocities, including 
spectroscopic binaries with centre-of-mass RVs, colour-coded the 
same as in Fig. \E\ The spherical polar axes have a radius of 900 
pc, centred on the Sun. The viewer's angle is the same as in 
Fig. |3] chosen to illustrate the 63° between the Earth's equator 
and the Galactic plane (longest horizontal line) , emphasizing that 
the COR AVEL giants are only visible from the Earth's Northern 
hemisphere. 

arcsec yr~^) and apparent magnitude [the most comprehen- 
sive photometric coverage of ob served RAVE sta rs is the 
2MASS Point Source Catalogue (jCutri et al.ll2003l ) so here 
we choose 2MASS Ks] in the equation 

Hks = ifs + 5 log ^ + 5 = Mks + 5 fog vt - 3.379, (13) 

which also expresses Hks in terms of the more fundamental 
properties of absolute magnitude (Mks) ^-nd tangential ve- 
locity [vt in km s^^). Fig. ll3l shows that the Hks~{J ~ Us) 
plane can divide the observed RAVE stars into three distinct 
stellar populations. 

The largest stellar population in the RAVE sample has 
J — Ks < 0.5 mag (~51 per cent). These stars are mainly 
thin disc F-G dwarfs [represented on the isochrone by the 
Main Sequence (MS) at Hks = 8-9 mag]. Although the sub- 
giant branch (represented on the isochrone by the horizontal 
branch along Hks ~ 6 mag) is a short phase of stellar evolu- 
tion (see the Hertzsprung Gap at J — Ks ~ 0.5 mag in Fig. 
I13p . there will be a non- negligible number of subgiants in 
the RAVE sample. Hks is unable to kinematically separate 
the subgiants from the F-G dwarfs because their values of 
absolute magnitude are too similar for their proper motions 
to resolve the overlapping kinematics. This results from the 
overlapping distance ranges of the two populations, due to 
the large spread in apparent magnitude in Fig. 1121 

The two populations with J — Ks > 0.5 mag have suf- 
ficiently different absolute magnitudes for their proper mo- 
tions to resolve the very different kinematics resulting from 
the very different distance ranges of the two populations. 
The fainter population (~5 per cent) are thin disc K-M 
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Figure 11. AitofT projection of the Galactic co-ordinates (in de- 
grees) of the 151 856 RVs in the internal RAVE data release (black 
dots) and the 24 748 publicly available RVs from the first data 
release (red dots). The vast majority of RAVE data is below the 
plotted celestial equatorial plane ((5 = 0°). The pattern is due to 
the survey field centres of the first and second input catalogues 
being defined on 5.7° and 5.0° grid spacing respectively (the field 
of view of the 6dF field plates is 5.7°). 
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Figure 13. Reduced proper motion diagram of 149 660 RAVE 
RVs (2196 either do not have 2MASS photometry or mea- 
sured proper motions ) . The solar-scaled Padova 2M ASS isochrone 
ijBonatto et al.l |2004| . red line) , chosen to be representative of 
thin disc stellar populations (2.5 Gyr old, initial metal fraction 
= 0.019, [Fe/H] Si —0.02 dex), has been given thin disc kine- 
matics (tiy fa 30 km s~^ in equation I13I I in order to interpret 
rave's stellar populations (approximately partitioned by the 
dashed lines). 
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Figu re 12. DEep Near-Infra Red Survey CDENIS. lEpchtein et all 
Il999l) 7-band selection functions of 131 632 RAVE RVs, ob- 
served with the two RAVE input c atalogues. Tycho-2 a nd the 
SuperCOSMOS Sky Survey (SSS, iHamblv et all |2001^ were 
used to generate the first and second RAVE data releases (see 
ISteinmetz et al1l2006l for details) and contribute 64 581 RVs to 
the plot (49 per cent). DENIS was used to define the current in- 
put catalogue, which is still being observed and will generate all 
future data releases and contributes 67 051 RVs to the plot (51 
per cent). (20 224 RAVE RVs are not included in the plot because 
of incomplete DENIS sky coverage.) 



dwarfs (represented on the isochrone by the MS at Hks 
— 9-14 mag). The brighter population (~44 per cent) are 
K-M giants (represented on the isochrone by the Red Giant 
Branch at Hks = 1-4 mag). Although the RAVE giant sam- 
ple is still dominated by the thin disc, their intrinsic bright- 
nesses probe distances that sample a statistically significant 
number of thick disc stars and a non-negligible number of 
inner halo stars. 



4.2 Distances to RAVE stars 

Fig. [13] shows that Hks bifurcation can be used to kinemat- 
ically select K-M dwarfs, without spectroscopically derived 
log g, to derive photometric distances. Their atmospheres 
are dominated by saturated molecular bands so their lu- 
minosity is not very sensitive to metallicity. This means 
isochrones for this population are very similar for different 
metallicities and ages, leading to small photometric distance 
errors. RAVE KO dwarfs probe ~50-250 pc from the Sun. 
This volume has already been sampled by the CORAVEL 
dwarfs (see Fig.O and shown to be devoid of tidal streams in 
Section (2] Therefore there is no advantage to searching this 
volume again with RAVE's less accurate W space velocities. 

The largest uncertainty in deriving photometric dis- 
tances to the RAVE F-G dwarfs is stellar age. Absolute mag- 
nitude is most sensitive to age close to the MS turnofT. Unlike 
the CORAVEL dwarfs, RAVE does not have an independent 
source of trigonometric parallax from which distances can 
be inferred. This enables stars to be placed in the H-R dia- 
gram and allows age to be estimated from a star's evolution 
away from the zero-age MS. In principle, a chromospheric 
activity indicator is present in RAVE spectra in the form of 
Ca II emission. This feature declines with time (as stellar 
rotation decreases with age) but decays into invisibility at 
about the age of the Sun. Ranking age using this diagnostic 
is fraught with systematic errors. RAVE F-G dwarfs probe 
similar distances (~200-400 pc) to the CORAVEL giants but 
are complementary in their non-overlapping sampled vol- 
umes (CORAVEL giants are Northern celestial hemisphere 
only and RAVE is Southern celestial hemisphere only - see 
Fig. [To}. Therefore, there is potential advantage to deriving 
RAVE F-G dwarf space velocities. However, the associated 
errors will not yield a reliable W distribution to search for 
asymmetries. 

Absolute magnitude is sensitive to [a/Fe] as well as 
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[Fe/H]. RAVE K-M giants include thick disc and halo stars, 
which have, compared to the thin disc, enhanced [a/Fe]. 
Metallicity uncertainty dominates the uncertainty in deriv- 
ing photometric distances to the RAVE K-M giants. 

Stars in the first RAVE data release were observed with- 
out a blue light blocking filter (OG531). Their spectra were 
contaminated with second order light, which means stellar 
parameters are not c urrently available for these stars (see 
ISteinmetz et al.ll200^ for details). Stars in all future RAVE 
data releases have all been observed with an OG531 filter 
so their spectra only contain first order light. Stellar param- 
eters can be derived for these spectra from the best-match 
synthetic spectrum. At the time of writing, the accuracy 
of these parameters is being finalized ready for the second 
data release (Zwitter et al. 2007, in preparation). The log g is 
likely to be sufficiently accurate to resolve the J — Ks > 0.5 
mag sample into dwarfs and giants. However, the log g er- 
rors may not be small enough to resolve the J — Ks < 0.5 
mag sample into dwarfs and subgiants. 

RAVE K-M giants probe greater distances (~0.7-3 kpc) 
than the CORAVEL giants. Therefore, there is great advan- 
tage to deriving the space velocities of RAVE K-M giants 
because it is unexplored phase space on Galactic scales in 
between the very local thin disc and inner halo at the so- 
lar position. However, tangential velocity errors scale with 
distance from the Sun. The absolute magnitudes of nearby 
Hipparcos red clump giants show no correlation with [Fe/H]. 
They have a mean Mk = —1.61±0.03 and ctmk ~ 0.22 mag 
l|Alvesll2000l ). which translates to a distance error of ~11 per 
cent. At 1 kpc, this distance error converts Afj, = 2 mas yr~^ 
to A«T = 10 km s"^ (Veltz et al. 2007) and ergo AW > 10 
km s~^. The red clump RAVE giants have the smallest dis- 
tance error of all the RAVE giants but their values of AW 
are still too large to resolve any dynamically cold, vertical 
tidal streams falling through the sample. The availability of 
RAVE giant metallicities would therefore not improve the 
situation. Outside the red clump, a RAVE K giant at a dis- 
tance of 1 kpc with proper motion errors at a more typical 
level of 3.5 mas yr~^ and a distance error >15 per cent, due 
to the uncertainty in [Fe/H], results in Avt > 20 km s~^ 
and ergo AW > 20 km s^^. 

4.3 RAVE radial velocity determination and 
accuracy 

The RAVE RV pipeline matches each RAVE spectrum with 
a theoretical spectr um from the IZwitter et al.l l|2004l ') and 
iMunari et al.l l|2005l ) libraries. These libraries contain thou- 
sands of stellar templates, whereas CORAVEL uses only one 
template. The RAVE RV pipeline uses the standard cross- 
correlation procedure l|Tonrv fc Davis! 19791), im plemented 
in the IRAF package XCSAO ("Kurtz et al]|l99d ). to derive 
RVs and internal RV errors. Fig. [14] shows the mode of the 
internal RV error distribution is ~1 km s~^. A RV zero-point 
error of ~1 km s~^ makes the mode of the external RV error 
distribution ^2 km s~^. 

We use the Tonr y-Davis cross-correlation coefficient (r, 
iTonrv fc Da^l 19791 ) to impose a quality cut. r is a measure 
of template match to spectra. Excluding RVs with r < 15, 
removes poor-quality spectra and a small number of hot 
dwarfs due to template-mismatches. This is seen in the pro- 
portion of RAVE RVs with J - Ks < 0.5 mag and r > 15 
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Figure 14. Distribution of internal RV errors for all 151 856 
RAVE RVs (all values of r, solid line) and the 140 228 RVs with 
r > 15 (dashed line). 



decreasing from ~51 to 50 per cent. Fig. 1141 shows this cut 
decreases the amplitude of the tail of the internal RV error 
distribution, which has sufficient accuracy to resolve dynam- 
ically cold streams. 

RAVE RV accuracy compares favourably with the 
CORAVEL accuracy because the two strongest Calcium 
triplet absorption lines in the RAVE spectroscopic region 
(8410-8795 A) are powerful RV diagnostics at RAVE's 
medium resolution {R ^ 7500). ISteinmetz et al.l (|2006l ) 
cross-matched RAVE first data release targets with the 
CORAVEL dwarfs. Thirteen matches were found but two 
of them were classified by N04 as binaries. The remaining 
eleven single targets show good agreement, with a mean dif- 
ference of 1.4 ± 0.4 km s~^ and a = 1.4 km s~^. 



4.4 RAVE pseudo-W without distances 

A star's heliocentric RV cannot be decomposed into the 
star's constituent U , V and W space velocities because, by 
definition, RV does not include any tangential contribution. 
However, a star's RV can be resolved into its line-of-sight 
(los) components of these space velocities: 



RV = ^(f/io<<)2 + {y^osy 4. {w^o,)2^ 



(14) 



We remove the dependence of Galactic latitude on RV by 
only resolving the components in equation [14] in a single 
dimension in the Galactic plane. 



RV cos 6, 



(15) 



and the dimension perpendicular to the Galactic plane, 

W'°"=RVsin6. (16) 

The above equations show that for a stellar los along b = 
—45°, RV consists of equal contributions from UV and 
W''°" but for stars with b < —45°, as b decreases the con- 
tribution of W'-"" to RV increases and RV becomes less and 
less sensitive to UV''°''. If a star's los is exactly towards the 
SGP (6 = -90°), W'^"" = W because here is there is no 
tangential component of W . For los with —90° < 6 < —45°, 
W^"" < W but W'"" is a useful proxy for W. Therefore, the 
issues discussed in the previous section can be circumvented 
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Figure 15. Lambert polar (equal area corresponds to equal solid angle on the sky) projections of the Southern Galactic hemisphere, where 
the SGP (b = —90°) is at the centre and concentric circles correspond to constant b at —60°, —30° and 0° (Galactic equator), around 
which the radial lines of constant I are labelled. CORAVEL dwarfs are plotted with their RVs in km s~^ colour-coded in different reference 
frames according to the legend of each plot: heliocentric RV (top left); peculiar RV of the Sun relative to the LSR (solar motion, RVpec, 
top right, see Appendix A for more details); RV with respect to the LSR (RV^sfl, middle loft); peculiar RV of the LSR relative to the GSR 
(RVpec, GSfli middle right); W''°" (RViSijsinft, bottom left). The overdensity at I ~ 180°, b ~ —20° is the Hyades open cluster. Bottom 
right: RAVE stars colour-coded according to RAVE field symmetry about I = 180° . The magenta curve traces the celestial equatorial plane 
(5 = 0°). 
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by using the accuracy of the RAVE RVs towards the SGP 
with h < -45°. 

The top two plots in Fig. [15] show that for differing 
los, RV in equation [16] is sensitive to differing amounts of 
solar motion with respect to the LSR. In order to compare 
along different los in the same reference frame, RV 
in equation [16] is replaced by RVlsa, which is RV corrected 
for the solar motion, decomposed into its cardinal directions 
along the los relative to the LSR: 

RVlsa = RV + U'^gn cos I cos h + sin / cos b 

+M/°s«sin6. (17) 

The middle left plot in Fig. [15] shows that the LSR reference 
frame removes the peculiar components of the solar RV (top 
right plot) from the RV los (top left plot). However, like RV, 
RVlsa is still dependent on the los. The middle right plot in 
Fig. [15] shows that RVlsb is sensitive to differing amounts 
of the peculiar RV of the LSR relative to the GSR, 

RVpec,GSH = Vrot siu / COS b. (18) 

The net Galactic rotation from / = 270° to I = 90° is just 
visible in the middle left plot of Fig. 1151 towards [l = 90°) 
and against (I = 270°) the Galactic rotation directions: 
the Q° < I < 180° CORAVEL dwarfs have sUghtly more 
stars with RVlsh < than RVlsh > and vice versa for 
180° < I < 360° (RViSij exhibits the pattern in equation 
1181 and the middle right plot in Fig. [T^ because it mirrors 
RVgsh). Although F05 showed that the presence of dynam- 
ical streams causes there to be a net radial motion in the 
solar neighbourhood (see Appendix A), this is in phase space 
and independent of I. Thus the 1/'°" contribution to HVlsr 
is isotropic as illustrated by the symmetry in RVlsa be- 
tween the regions near the cardinal U directions {I — 0° and 
I = 180°) in the middle left plot in Fig. [15] 

It is not appropriate to correct WJlsr for Galactic ro- 
tation because we are resolving RVlsh perpendicular to the 
Galactic plane and the correction does not affect W''°" be- 
cause W — in equation 1181 However, because this cor- 
rection is not applied to RAVE stars towards the SGP 
(fe < —45°), the small component of RVlsb sensitive to 
Galactic rotation (1/'°") propagates through equation 1161 
The effect of V''"" is just visible in the bottom left plot 
in Fig. 1151 where sin(6 < 0°) reverses the sign of the ef- 
fect i.e. the 0° < / < 180° CORAVEL dwarfs have slightly 
more stars with > than < and vice versa 

for 180° < I < 360°. A Kuiper symmetry test on the W'°° 
distribution of & < —45° CORAVEL dwarfs is similar to the 
results in Section 12.41 because the W''"" sample is symmet- 
ric about I = 180°, so that the positive and negative 1/'°" 
contributions cancel out. 

A Kuiper symmetry test on the W^g^ distribution of 
RAVE stars with b < —45° and r > 15 strongly rejects the 
null hypothesis (that Wtin > and W^in < are drawn 
from the same parent population) at Q ~ 1 x 10~®, where 
W'j°g'^ = 4.1 km s~^ to the nearest 0.1 km s~^. The bottom 
right plot in Fig. 1151 shows that, because the RAVE sample 
is Southern celestial hemisphere only, it does not symmet- 
rically sample the SGP out to 6 = —45°. This plot reveals 
the largest part of the b < —45° asymmetry at I > 180° 
(colour-coded green) is not observed by RAVE at / < 180° 
because it is in the Northern hemisphere. The contribution 



of 1/'°" to KVlsr in this region of I > 180° cannot cancel 
with its I < 180° mirrored region so the Kuiper symmetry 
test finds the difference in the Wj^gji distribution. 

To rectify this, we only sample b < —45° RAVE fields 
that are symmetrically positioned either side of / = 180° 
(28 868 RVs colour-coded red in the bottom right plot of 
Fig. ll5|l . Because RAVE observes some fields more than once 
(normally different stars each visit) , differing stellar densities 
are visible in this plot. In the Galactic sky-symmetric sam- 
ple, there are more r > 15 RVs with I > 180° (16 235) than 
I < 180° (12 633). Nevertheless, the next section shows that 
the Kuiper symmetry test fails to reject the null hypothe- 
sis, showing that our Galactic sky-symmetric sample can be 
used to search for streams without any further processing. 
Its selection function is almost identical to Fig. 1121 apart 
from nearly all the 9 < 7 < 10 mag stars come from the first 
input catalogue rather than the second. The stellar popula- 
tion proportions of the subsample are different to the whole 
RAVE sample. F-G dwarfs have increased most to ~61 per 
cent, K-M dwarfs have increased to ~7 per cent and K-M 
giants have decreased to ~32 per cent. The internal RV er- 
ror distribution of the subsample is the same as the whole 
sample with r > 15. 

4.5 Determining RAVE W^sr 

We repeated the Kuiper test in Section [2.31 comparing the 
b < -45° RViSflsinb distribution of all the CORAVEL 
dwarfs with the single dwarfs and found them to be more 
similar than the Wlsr test (Q = 0.9998 but Ne is only 779). 
This demonstrates that all the RAVE (and SEGUE) RVs can 
also be used to look for tidal streams without worrying about 
the affects of binarity. This section repeats the technique 
applied to the CORAVEL stars in Sections [2]4] and [3]3] to 
RAVE stars. 

The top of Fig. [TUl shows the Wj^gji distribution for the 
minimum value of D, where W'j°g'^ = 1.7 km s~^ to the 
nearest 0.1 km s^^. The difference between this value and 
the W®gj^ = 7.0 km s"^ values derived from the CORAVEL 
dwarfs and giants highlights that Wj^gji ^ Wlsr but merely 
a proxy. The D- and D+ positions in the bottom plot of 
Fig. [T5] are similar to the CORAVEL positions in Figs [3] 
and [7] suggesting W[°gji does approximate Wlsr- W^j°g^ is 
nearly in agreement with the mode of the Wlsr distribution 
in Fig. 1161 Again, this suggests that there cannot be many 
stars belonging to vertically coherent streams with systemic 
Wlsr velocities. Table 5 shows that the RAVE A^e value 
is more than double the value for the CORAVEL samples. 
The RAVE Q value is double the value for the CORAVEL 
dwarfs but less than the value for the CORAVEL giants. 
The RAVE AW^^^jf values are more symmetric than the 
CORAVEL values. 

4.6 Sensitivity of RAVE W'lsr symmetry to tidal 
streams 

To estimate the volume probed by our RAVE sample, we 
consider that RAVE stars fainter than 7 = 12 mag contribute 
a negligible number of stars to the sample (see Fig. ll2p . This 
sets a distance limit for a typical RAVE giant {Mj — —0.25 
mag) at ~3 kpc. We approximate the volume of each RAVE 
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Table 5. Results of the Kuiper symmetry test, applied to our RAVE sample, to determine W^g'^ to the nearest 0.1 km s ^, 
where /^W]°^'® = W''°^^^ - 1.7 km s"!. 
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D+ 
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0.4 


-1.3 


14 742 


14 126 
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0.0318 
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0.0323 


1 X 10"^ 
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0.0146 


0.0053 


0.0199 


0.0658 


2.7 


1.0 


13 455 


15 413 


7184 


0.0119 


0.0196 


0.0315 


3 X 10-5 



Table 6. Results of the Kuiper symmetry test, applied to our RAVE sample, to determine the approximate number of tidal stream 
stars (Ns), required in pseudo-randomly generated <Jb = 10 km s"^ Gaussians, placed at ±(1, 2, 3, ~7)o'j^!os in the Wj^gj^ 

distribution, to cause the test to reject the null hypothesis at 4(t (Q < 6 X 10~^). % = Ng/{Ng + N-^^^^i) where N^^^^i = 28 868. pg 
is the range of stellar density of each stream, in N stars kpc" calculated by dividing the Ns of each stream by our RAVE sample 
volume (~8 kpc" ^, see text for details), modulated by the 5-15 per cent RAVE completeness, which increases Pa by a factor of 
20-6.7 respectively. 
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field as a cone of height 3 kpc subtending tiie RAVE field 
diameter of 5.7°. Our sample is complicated by the inclusion 
of stars from two input catalogues witii the co-ordinates of 
its field centres based on grids offset from each other. Fig. 
[IT] shows tiiat the first RAVE input catalogue (colour-coded 
red) consists of contiguous fields on a 5.7° grid. The second 
RAVE input catalogue (colour-coded black) is based on a 
5.0° grid and so the 5.7° field of view 6dF fields overlap 
each other, as well as overlapping the fields of the first input 
catalogue. We have counted the number of contiguous (red) 
fields in Fig. [17] and estimated, to the nearest half field, 
the additional sky covered by the non-overlapping parts of 
the other fields (black). The sample consists of ~113 fields, 
giving a volume of ~8 kpc'^. 

^Steinmetz et al.l l|2006l ) estimated the RAVE complete- 
ness of the old Tycho-2 and SSS input catalogue in their fig. 
4 to be ~15 per cent compared to DENIS at the bright end 
of the selection function (7 ~ 10 mag). This drops to ~5 
per cent at the faint end (J ~ 12 mag). Subtle colour biases 
exist in the old input catalogue but more than half our sam- 
ple comes from the unbiased DENIS input catalogue. This 
should reduce the effect of the biases and increase the com- 
pleteness levels but we use the old input catalogue levels as 
a conservative estimate. 

This section repeats the technique applied to the 
CORAVEL stars in Sections [23] and [Tito the RAVE stars, 
except the RAVE dispersion is much larger: f^^io^ « 29 km 

s"^ and thus now zbScr^yios approximately corresponds to 
the ±W]^sR values of the II99b streams. 

Table 6 shows that the positions of D- and D+ appear 
to cause the +W\^gji distribution to be insensitive to the 
position of the stream and the number of its members. As 
before, the test becomes more sensitive (less stream stars 



required to generate a Aa detection), the further a stream 
is from the centre of the distribution in the —Wj^gji distri- 
bution. However, Table 6 shows that, again, the number of 
stream stars required for a 4cr detection stays approximately 
constant (saturates) at Wlsr ^ ~i<^WLSR- 



Fig.[T5]shows that there are a few RAVE stars with sim- 
ilar Wlsr as both the +Wlsr (250 < Wlsr < 350 km s^^) 
and -Wlsr (-250 < Wlsr < -200 km s"^) H99b streams. 
Without distances to these stars to check their orbital an- 
gular momenta, they cannot be confirmed as members of 
the H99b streams. Regardless of the symmetry between the 
two streams. Table 6 shows that there are far too few of 
these stars for the Kuiper symmetry test to detect them as 
a stream. 



Assumi ng 5 per cen t RAV E completeness, the lower 
bound of the lrVeese et ahl (|2005l ) estimates of the Sgr stream 
stellar density predict that ~80 RAVE stars could be Sgr 
stream stars. This increases to ~590 stars using the upper 
bound. Assuming a less conservative RAVE completeness 
level of 15 per cent, predicts between ~250-1780 Sgr stars. 
The number of VOD stars that could be in RAVE is >48 
000. Therefore, the superior volume of RAVE rules out the 
possibilities that the Sgr stream and VOD cross the Galac- 
tic disc through the solar neighbourhood. It also puts much 
tighter constraints than the CORAVEL surveys on the den- 
sity of any stream that could be passing through the solar 
neighbourhood, ruling out any streams with similar densities 
to the Sgr stream. 



Searching for streams in CORAVEL and RAVE 19 




300 



270 



240 



180 



Figure 17. Same as the bottom right plot of Fig. [15] except only the Galactic sky-symmetric RAVE sample is plotted and is colour-coded 
according to input catalogue: 5.7° grid spacing (red) and 5.0° grid spacing (black). The blank strips (aligned with the declination axis) in 
between the black dots are due to the incomplete sky coverage of DENIS data. 



4.7 



Comparing RAVE and CORAVEL W^^'°|j' 
sampled volumes 



and 



Table 7 presents the Kuiper test results from comparing the 
entire RAVE W\^sr distribution with the entire CORAVEL 
dwarfs (Fig. [T9|l and giants Wlsr distributions (Fig. [20|) . 
The test finds the largest differences at ~±lcrwj,sji in 
the CORAVEL distributions because of the greater RAVE 
(j^yio^ , due to the larger volume it samples. A comparison 
of the respective volumes sampled by the surveys is depicted 
in Fig. [21] It shows that only the ~2000 RAVE K-M dwarfs 
probe the same volume as the CORAVEL F-G dwarfs. The 
~17 600 RAVE F-G dwarfs fill in the SGP at 6 < -45° 
not probed by the CORAVEL giants. It is the ~9200 RAVE 



giants that really extends our tidal stream search into the 
thick disc, its scaleheight being ~1 kpc (Veltz et al. 2007). 
Because of the differing distance ranges covered by the differ- 
ent RAVE stellar populations, we have repeated our search 
procedure on each population separately (as defined in Fig. 
I13p . This has been investigated in case e.g. the asymme- 
try of a thick disc or halo stream of giants is smeared out 
by the more numerous sample of dwarfs. However, a Kuiper 
symmetry test on each component fails to reject the null hy- 
pothesis. F-G dwarfs are found to be the most symmetric, 
then K-M dwarfs with the K-M giants the least symmetric. 
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Figure 16. Top: W|°^jj distribution of RAVE stars with h < 
—45° (solid histogram) either side of the Wj^gj^ = dividing line 
(dashed lines). The 5 km s~^ velocity bin sizes are chosen to be 
more than twice the mode of the RAVE external RV error (~2 
km s-l). Bottom: Normalized CDF of < (solid lines) 

l'sr ^ ^ (dotted lines), showing D- (dashed lines) and 
(dot-dashed lines). 



and 

D. 



Table 7. Kuiper test results comparing the RAVE Wj^^j distri- 
bution to the CORAVEL dwarfs and giants Wlsr distributions. 



CORAVEL 




D- 


D+ 


D 


Q 


Dwarfs 
Giants 


9077 
4988 


0.1060 
0.0863 


0.0974 
0.0538 


0.2034 
0.1401 


<io-»7 

<10-3^ 



5 DISCUSSION 

We have searched the CORAVEL and RAVE surveys for evi- 
dence of a net vertical flow through the solar neighbourhood 
that could be associated with a tidal stream of stars with 
vertically coherent kinematics. Due to the current unavail- 
ability of accurate distances to RAVE stars, we have used 
the Kuiper statistic to test the symmetry of the vertical ve- 
locity distribution function of the samples. In particular, the 
test serves as a useful first search for the presence of promi- 
nent, kinematically coherent, vertical tidal streams in the 
RAVE data, which probes further into thick disc and inner 
halo phase space than any other survey of the local Galactic 
volume. Using the Geneva-Copenhagen stellar binarity flag, 



Figure 18. Top: RAVE Wj°^j^ distribution (solid histogram) 
either side of the Wj^gj^ = dividing line (dashed line). The 
randomly-generated Gaussian (representing a cr = 10 km s~^ 
tidal stream) included in the data set (dotted histogram) that 
cause the Kuiper symmetry test to reject the null hypothesis at 



4(T when placed at ~— So-pj/^g^ is also plotted separately (dot- 
dashed line). Bottom: Normalized CDF of VK|°^=^ < (solid line) 
and Wj^gj^ > (dotted line), showing D_ (dashed line) and 
(dot-dashed line, just visible at lVK|^'g^| ~ 0). 



we demonstrate that binarity has a statistically insignificant 
effect on the W and RViSflsinfe velocity distribution func- 
tions of the CORAVEL dwarfs, allowing us to use all the 
CORAVEL dwarfs and RAVE radial velocities (RVs) in our 
search. 

A tidal stream falling slowly through the solar neigh- 
bourhood volumes considered (see Table O needs to be very 
coherent for our Kuiper symmetry test to detect it against 
the background W^gn^ distribution at |W];s|;'| < o- iias). 

As expected, there is a general trend that the test can de- 
tect smaller numbers of stars in a coherent stream the fur- 
ther from the centre of the distribution the stream is, out 



to ~±3ct. 



However, further into the tails of the dis- 



tribution the test becomes insensitive to the position of the 
stream and it is the number of stream members that cause 
the and — W'^glj^ CDFs to diverge sufficiently to 

generate a stream detection. This number appears to ap- 
proximately saturate at |Vl^^g|;'| > 3(j^,(io3) and is set by the 

detection threshold at ~±3a 



^ LSR 



LSR 

) . The required number of 
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Figure 19. Top: W^gjj' distribution of the 28 868 RAVE stars 
(solid histogram) and the 13 240 CORAVEL dwarfs (dotted his- 
togram). Bottom: Normalized CDF of RAVE (solid line) and 
CORAVEL dwarfs (dotted line) as a function of W^g^'j^, where 
the maximum differences between them are indicated by vertical 
lines: (dot-dashed line) and D_ (dashed line). 
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Figure 20. Top: VK^gJj' distribution of the 28 868 RAVE stars 
(solid histogram) and the 6030 CORAVEL giants (dotted his- 
togram). Bottom: Normalized CDF of RAVE (solid line) and 
CORAVEL giants (dotted line) as a function of VF^'g^j', where 
the maximum differences between them are indicated by vertical 
lines: Dj^ (dot-dashed line) and D_ (dashed line). 



stars is >200. Fig. [18] shows streams should be detectable 
by eye in the extreme tails of a distribution where more 
specialized statistical analysis could be applied to test their 
significance. We have not found any asymmetries in any of 
our samples that warrant such treatment. Table 8 summa- 
rizes the limits we have been able to place on the properties 
of streams in the solar neighbourhood. 

I n contrast to the 'Fiel d of Streams' found in the outer 
halo l|Belokurov et al.]|2006l ) , we find that the volume of the 
solar neighbourhood sampled by the CORAVEL and RAVE 
surveys is devoid of any tidal streams with coherent verti- 
cal kinematics containing hundreds of stars. This suggests 
the Sagitarrius (Sgr) tidal stream does not pass through the 
solar neighbourhood, whic h agrees with the late st empirical 
determination of i ts orbit (jNewberg et al.|l2007l) and disso- 
ciates it from the iHelmi et al. stream l|Freese et al.l 

|2004| ). 

The absence of the Sgr stream near the Sun is consistent 
with simulations of the disruption of Sgr in near l y spherical 



and pr olate Galactic p ot ential s (|Helmil 20041: 
iFellhauer et al.l 



2005 



Law et al 



l2006l : iMartmez-Delgado et al 



20071 ) and seemingly inconsistent with oblate potentials 



(|Law et all l2005l : iMartfnez-Delgado et aLll2007l ). The con- 
straints on the number of the Sgr stream stars in the solar 
neighbourhood volume searched in our samples could prove 
a useful tool for discriminating between Galactic potential 
models. 

The CORAVEL Kuiper symmetry test results are sup- 
ported by the lack of strong clumping in orbital angular 
momenta space. The limited range of the CORAVEL g iants 
in angular momen ta space shows the iMaiewskil (|l992l ) and 
iHelmi et al.l (jl999l ). hereafter H99b, streams are not present 
in this sample. One of the CORAVEL dwarfs is a member of 
one of the H99b streams. There is another CORAVEL dwarf 
with angular momenta and metallicity consistent with the 
H99b streams. However, we consider it to be a possible out- 
lier because its kinetic energy is too large to be consistent 
with the energies of the other members of the H99b streams. 

The lack of a net vertical fiow through the solar neigh- 
bourhood argues against the Virgo overde nsity (VOD) cross- 
ing the disc near the Sun. This agrees with ljuric et al.l (|2005| ) 
not finding it within a few kpc of the Sun. Their preliminary 
analysis of 2MASS M giants did not reveal a similarly large 
density enhancement to the VOD in the Southern Galactic 
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Figure 21. Three dimensional Galactic sky distribution of all 
the 13 092 CORAVEL giants with W space velocities (grey), the 
6030 CORAVEL giants with W space velocities, including spec- 
troscopic binaries with centre-of-mass RVs, colour-coded accord- 
ing to their maximum likelihood base group assignment in Fig. [B] 
and the south Galactic cap RAVE sample (orange), where each 
RAVE star has been assigned a distance of 300 pc, which par- 
titions the maximum distance range of RAVE K-M dwarfs and 
the minimum distance range of RAVE F-G dwarfs, and a dis- 
tance of 3 kpc to represent the maximum distance range of / = 
12 mag RAVE K-M giants. The spherical polar axes have a ra- 
dius of 3.7 kpc, centred on the Sun, chosen to include the most 
distant CORAVEL giant with a W velocity (3.6 kpc, marked by 
the black cursor) with the same viewer's angle as in Fig. 1101 



Table 8. Summary of how many stream stars {Ns) and their den- 
sity (ps in N stars kpc~^) are needed for the Kuiper symmetry 
test to detect them and how many stars (A^) from the Sagitarrius 
stream (Sgr) and Virgo overdensity (VOD) are expected in our 
sample volumes (V), where the CORAVEL values correspond to 
the volumes within which they are complete and the RAVE vol- 
ume is its total approximate volume (its volume completeness 
(VC) is taken into account). Whether the Kuiper symmetry test 
is sensitive to the presence of each stream, given their non-local 
stellar density estimates in the literature is indicated: yes (y), no 
(n), maybe (?). 



Sample 


CORAVEL 
Dwarfs Giants 


RAVE 




Section 


[231 


13:41 


\m 




V (kpc3) 


0.0003 


0.0511 


7.9052 




Ns (low) 


200 


200 


300 




Ns (high) 


600 


800 


600 




VC (%) 


100 


100 


5 


15 


ps (low) 


0.7 X 10*^ 


4000 


800 


300 


ps (high) 


2.2 X 10'^ 


16 000 


1500 


500 


N Sgr (low) 


0.1 (n) 


10 (n) 


80 (n) 


250 (?) 


N Sgr (high) 


0.4 (n) 


80 (n) 


590 (y) 


1800 (y) 


N VOD 


30 (n) 


6000 (y) 


48 000 (y) 144 000 (y) 



Since Jz is conserved duri ng phase-mixin R , tida l streams 
should have similar Vgsr- iGilmore et al.l (|2002l ) found a 
tidal stream with coherent Vgsr in between canonical thick 
disc and inner halo values. The magnitude range of RAVE 
was chosen so that its giants probe the interface of the thin 
and thick discs and inner halo. Hence, RAVE giants can be 
also used to investigate the angular momenta of the Milky 
Way's stellar components and its intervening stream(s?) 
(Seabroke et al. 2007, in preparation). 



hemisphere. All this evidence is in accord with dark matter 
being smoothly distributed in configuration space through- 
out the solar neighbourhood. 

Comparison between the vol ume covered by SP SS pho- 
tometry in the top right plot of lJuric et al.l (|2005l ) fig. 20 
and the local volume with measured pha se space in Fig. 121! 
highlights the disparity between the two. ljuric et al.l (|2005l ) 
extrapolate the four thin and thick disc overdensities seen 
in their survey volume to the full Galactic disc {\Z\ < 3 kpc, 
7? < 15 kpc) to imply that there are ~20-40 substructures of 
this type in the Galaxy, which could be the "missing satel- 
lites". It is not too surprising that none of these are within 
the small Galactic volumes searched in this paper. 

Depending on the availability and accuracy of metallici- 
ties in the second RAVE data release (Zwitter et al. 2007, in 
preparation), distances to RAVE giants may be sufficiently 
accurate to allow action-angle space or integrals of motion 
space to be searched in the future for clustering indicative 
of tidal streams. Addi tional photometric data in the future 
(e.g. from SkvMapper. lKeller et al]|2007h will further reduce 
distance errors to RAVE stars. However, as in this study, ac- 
curate metallicity-derived distances are not required if using 
RAVE stars towards the cardinal Galactic directions. In the 
rotation direction (Z ~ 270°), RVs are dominated by Vgsr, 
where the conversion between the t wo is only weakly de - 
pendent on distance (|Woollevlll978l : iFrenk fc Whit3ll980^ . 
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APPENDIX A: HISTORICAL PARALLELS 

Soon after IVogej (ll873l'FI measured the first reliable RVs of 



Sirus and Procvon. ISeabrokd (|l879l fl performed one of the 
first RV surveys, measuri ng 68 RVs for 2 9 stars, followed 
by 699 RVs for 40 stars (ISeabrokd [l887al lbh and 866 RVs 
for 49 stars (|Seabrokel Il889l ). Although these visually de- 
termined RVs were individually inaccurate, when analysed 
statistically they showed that the Su n was moving r elative to 
nearby stars at 24.5 ± 15.8 km s"-' l|Homannlll88t j'). Within 
the (admittedly large) error, this measurement of solar speed 
with respect to the LSR agrees with the standard modern 
measurement of 13.38 ± 0.12 km s~^ (calculated from the 
iDehnen & BinnevI (119981 ) solar sp ace velocities in Table Al). 
Fig. lAll shows the direction of the (jHomann!ll886l) solar apex 
towards I = 4 4.1°, b = 8.5° is approxim ately in the same di- 
rection as the iDehnen & BinnevI (|l998l ') values: I = 27.72°, 
6 = 32.43°. 

We have resolved the iHomannI (|l886l ) solar motion into 
its Galactic cardinal direction vectors, t^®g^, and 
W^gj^, assuming the solar apex direction is error-free. It 
almost certainly is not but iHomannI l| 18861 ) does not quote 
a directional error so the solar space velocity e rrors i n Table 
Al should be considered a lower limit. The Homann (1881) 



U 







LSR agrees with more modern estimates within the errors 
while V^gj^ and W^^j^ are close to agreeing within the er- 
rors. 



Over a century after the iHomam] (|l886l l derivation of 
the solar motion, the debate on this subject in the litera- 
ture has been revived by the CORAVEL surveys. F05 in- 
terpreted the solar neighbourhood substructure in the U- 
V plane found in their sample of CORAVEL giants and 
in CO RAVEL dwarfs by N 04 (and in Hipparcos tt, jj, and 
RVs bv lSkulian et al.lll999l ) as dynamical streams. Table Al 
shows that when F05 excluded all the dynamical streams, 
they found a different C^fg^ from when t hey included all the 
CORAVEL giants, which agrees with the IDehnen fc BinnevI 
1 19981 ) value. This discrepancy raises the question of how to 



^ A predecessor as director of Potsdam Observatory to MS, co- 
author of this paper and PI of the RAVE project. 

The great-great-grandfather of GMS, first author of this paper. 
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Table Al. Comparison of the iHomannI l|l886h solar space velocity components compared to more modern 
estimates. 



Author(s) 


Data used 


^ LSR 
(km s^-"-) 


^LSR 
(km s-i) 


WfsR 
(km s ^) 


Homann fl886) 


Seabroke (1879) RVs 


17.4 ± 11.2 


16.9 ± 10.9 


3.6 ± 2.3 


Dehnen & Binnev ( 1998) 


Hipparcos it &i fi 


10.00 ± 0.36 


5.25 ± 0.62 


7.17 ± 0.38 


This paper 


CORAVEL dwarfs 






7.0 


This paper 


CORAVEL giants 






7.0 


F05 


CORAVEL giants 


10.25 ± 0.15 




7.98 ± 0.09 




Excluding streams 


2.78 ± 1.07 




8.26 ± 0.38 


Veltz et al. (2007) 


RAVE RVs, UCAC2^ ^i 


8.5 ± 0.3 




11.1 ± 1.0 



" UCAC2 is the Second US Naval Observatory CCD Astrograph Catalogue l|Zacharias et ahlliooi ) 
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Figure Al. Lambert polar (equal area corresponds to equal solid 
angle on the sky) projections of the entire Galactic sky, where the 
SGP (6 = —90°) is at the centre and concentric circles correspond 
to constant b at every 30° to the outermost circle at +90° (NGP), 
around which the radial lines of constant / are labelled. The Galac- 
tic directions of estimated values of the solar apex (triangles) and 
antapex (squares) are plotted. The thick da shed line connects 
the so lar motion direction calculated from the iDehnen &: BinnevI 
lll998l ) solar space velocities in Table Al. The thin dashed lin e 
connects the solar motion direction calculated bv lHomannllll886h . 
calculated from lSeabrokS l|l879h RVs. 



motion, unbiased by the presence of in-falling tidal streams 
on to the Milky Way disc near the Sun. 

This paper has been typeset from a Ti^X/ JM^jiK file prepared 
by the author. 



derive the solar motion in the presence of dynamical pertur- 
bations altering the kinematics of the solar neighbourhood. 
They posed the question of whether a subset of stars ex- 
ists in the solar neighbourhood that has no net radial mo- 
tion, which can used as a reference again st which to mea- 
sure t he solar motion. The simulations of iDe Simone et al.l 
demonstrate that transient spiral waves can excite a 
net radial motion through the solar neighbourhood of ~10 
km s~^. 

Table Al also shows that F05 found that the presence 
of dynamical streams in the solar neighbourhood has a neg- 
ligible affect on W^^^. This suggests that the dynamical 
streams do not have a net vertical motion because they are 
well phase-mixed. In this paper, we have shown that there is 
no net vertical flow through the solar neighbourhood, sug- 
gesting is an appropriate measure of the vertical solar 



